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ABSTRACT 


An analytical study of stresses in the blades recorded during the March 1971 
tests of the DH 2011 jet flap rotor has been performed and presented. The main 
objective of the study was to compare the experimental results with analytlcaly 
determined stresses. The comparison extended over 15 specific flight cases has 
been only partially successful. In fact computed 3P and 4P stress components 
showed only a poor correlation with the test data obtained. It is believed that the 
simplified model of aeroelastic effects used is mainly responsible for this 
lack of agreement with test results. 
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1. INTRODUCTION 

The March 1971 Wind tunnel tests of the DH 2011 jet flap rotor have showed 
high correlation between stresses and multicyclic control effects. By a suitable 
choice of these multicyclic effects, reductions of 50% in stresses and vibrations 
have been obtained. Though these results correspond to what could have been 
expected, no reliable analytical methods exist to predict and explain the expe- 
rimental phenomena. The purpose of the present study is to attempt with existing 
analytical means the determination of stresses and their comparison with recor- 
ded data. 

The existing means consist of the computer program (Ref. 1), Icnown as the 
Evans - Me Cloud Program, EMCP , and an anal 3 dical method for determination 
of stresses and vibrations on the rotor blades from known aerod 3 mamic load 
distributions (Ref. 2). The Glravions Dorand Stress Program, GDSP, is based on 
transfer functions and transfer matrices suitable for investigation of complex 
rotary wing problems. 

The stresses will be analyzed at selected flight configurations, mainly at 
advance ratios of 0.4, which have been used in the stress and vibration reduc- 
tion study (Ref. 3). 


2. SYMBOLS 


Symbols of Ref. 1 are used throughout this document and will not be restated, 
•j" period 

complex Fourier coefficient of n-th order 
OC ^/t non dimensional time 

N number cf blade elements 

tip mass of n-th blade element 

/ 

m ^ central mass cf n-th blade element 

equivalent spring effect between n-th and n + 1st blade element 
i ^ length of n -th blade element 


angle between n-th and n + 1st blade element 
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shear at the end of the n-th blade element 
C bending stress 

aerodynamic force applied on m fi, 

^ ordinate at the end of the n th blade element 

angle between the n th blade element and the reference line 
E I local rigidity of the blade 

Am. mass distribution of the blade (blade density) 

dr 

I blade longeron thickness 

R rotor radius 

J1 rotation speed (rd/s). 


Upperscripts ( x ) and ( x ) indicates first and second order derivatives 


and 


* 


d X 
d t 


3. BASIC TEST DATA 

The basic stress data has been processed to obtain their analytical content. 
The followir^ method has been adopted ; the stress signals have been originally 
recorded on photographic strips of CEC recorders (fig. 1), The period then has 
been divided into 50 to 60 intervals defined by X , the corresponding ordinate 
y was noted. Traces n®*s 19 and 20 corresponding respectively to the stress 
at 45% radius (4. 2)and 70% radius (4. 3), have thus been treated for 73 test points. 

The Fourier analysis of the stress signals has been obtained at the Ames 
Research Center by a standard computer routine. 

In fact, only 15 c£ these analyses are of use for this study, since we have only 
15 simulations. The simulations are defined such that the rotor is in the same 
aerodynamic configuration as in certain windtunnel tests. 

In particular, the cyclic and multlcycllc pitch law was defined by Fourier 
analysis cf the flap deflection signal for these same tests. The deflection, produ- 
ced by a CIMATRAN transducer, is given by trace n® 24 of the C.E. C. recordings. 


O' 
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These signals were first analyzed by Giravions Dorand, using a mechanical 
harmonic analyzer and were then reanalyzed more accurately by NASA, using a 
computer. 

We therefore possess 15 simulations giving the radial distribution of forces on 
the blade as a function of the azimuthal position of the blade. These simulations, 
numbered SIM 11 through SIM 25, correspond respectively to run-points 9-3, 9-4, 
9-5, 9-6, 12-10, 12-11, 12-12, 12-13, 14-10, 14-11, 14-12, 14-13, 16-8, 16-9 
and 16-10, for which we know the Fourier coefficients of stress at 0. 45Rand 0. 7R. 


4. CALCULATION PROCEDURE 

The program described in the present note has been motivated by the need to 
correlate the theoretical results obtained from the Evans - Me Cloud simulation 
program (Ref. 1) with the experimental data gathered during the March 1971 wind 
tunnel tests of the DORAND jet flap rotor. In the E. M. C. program, the aerod 3 ma- 
mic forces are computed at several stations on the blade and for several blade 
azimuthal angles V . On other hand, blade stresses have been measured at 0. 45 
and 0. 7R, as well as the vibratory forces on the rotor hub, as a function of 
(Ref. 2 and 3). 

Thej purpose of the Giravions Dorand Stress Program is to compute the blade 
stresses and vibratory forces from the aerodynamic forces. The method is based 
upon the computation of transfer functions obtained from the dynamic and elastic 
characteristics of the blade (Ref. 4 and 5), which are then used to compute the 
Fourier coefficients of the stresses from the Fourier coefficients of the forces. 

The G. D. Stress Program is composed of three independent programs. 

1. BFHA (Blade forces harmonic analysis) 

Computes the aerod 3 mamic forces Fourier coefficients from the values of 
these forces given for discrete values of 'V 

2. BSTM (Blade stresses transfer matrices) 

Computes transfer matrices parameters from blade characteristics. 

3. STRESS 

Computes the Fourier coefficients of blade deflections and stresses and 
of the vertical vibratory force on the rotor hub, using the output cards of 
BFHA and BSTM. 

In this chapter of the present note the mathematical methods involved in these 
programs will be described. The following chapter will be devoted to the description 
and use of the FORTRAN IV programs. 


pNQ-t OG.A V ■ IVQp si 
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4.1 The fast Fourier analysis {Program FHAKM) 
4.1.1 The complex Fourier serie 


A periodic function F ( t ), of period X , may be expanded in a com- 
plex Fourier serie 


F (t) 



Zrr\ nfc /X 

a 


(1) 


The complex coefficients C are given by ' 
i 

- 2 TTJ nsc 

f (cc) t dac. 


where F {Toe) 

and oc = ^ / T 


(2) 


These coefficients are related to the coefficients of the cosine/sine 
expansion by the relations : 


r '’o 


= 2 X Rtt-oL (tort (C^) 


(3) 
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The computation of the Fourier coefficients amounts therefore to 
that of the integrand in (2), which may be approximated, in the case 
of digital computation, by the sum : 

z: (4) 

P=1 

The value of F is assumed to be known for N ^ points of the in- 
terval ( 0 , X ), so that 

f ^ ) = F Ct p ) 

O 

^ 1 


Fp = 

oc, = 

^ H+1 
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4.1.2 Principles of fast analysis 

A direct computation of the sum in (4) is not desirable because : 

1) it replaces the function f by a series of steps 

2) the exponential FORTRAN IV function has to be called up for 
each p and each r\ , which is a time consiuning process. 

The method of the fast analysis will consist in replacing the actual 
function f by linearly interpolated segments between two known 
values Fp and F p+i , and to compute the exact solution for the 
Fourier coefficients of this interpolated function f * 

Also, a recursive algorithm is used which deduces the value of 
the exponential function needed at a particular step from its value 
at the previous step by simple multiplication, using the fact that ; 

The FORTRAN IV exponential function is therefore called up only 
once, at the beginning of the program, and subsequent calls are 
replaced by a multiplication, which Is some 10 times faster. 

4. 1. 3 Equation of the fast analysis 


The function F is presumed known at equidistant points OCp . If 
D is the interval between two points : 






SCp . ( p-1) JD (6) 

Equation (2) may be written as : 
n'' 

” r _ s tt j n. ac 

* X / F C«) a- dae (7) 

pel •^pml 

Between aCpand ■Xp.'i , the function f is approximated by the linear 
function : 


f Ccc) = f, 


Fp+’ - fi 


- OCp 


(oc,ocp; 
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Replaclng f by p in Eq. 7, the integrand can be evaluated analy- 
tically and Cn is found to be equal to : 


_1_ [A 

2TTjn L 




where 


£rv = 


-"a- 


8n = 


_ 2 TT j a D 
a 

2irjn.D 

) / ZnjnD 




Since ft D = 1, it follows that 


=f, 


Therefore 


P-1 

fp = f, 


^ P+l ^N' 


Hence (9) may be written as : 

” 2^ ^ B„./ e„.j "Z^ fp 

or, using definitions (10) : 


1 — Park C 

-rr- 2- 

TT D 


^ fp£rv"'' 
p=T iP ^ 


This is the formula which is used in the program for computatii]^ 
the Fourier coefficients. 

The quantities £ ^ are obtained recursively by : 

Eri= Ei ^n.-^ (12) 
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The sum over p , by the recursive formula 


3 ^ = O 


giving 




p_i 


Finally, for n=o , a separate computation is performed : 


, 2 The method of transfer 
4.2.1 Principle of the method 

In this method, the blade is replaced by an equivalent system of N 
rigid elements linked to each other by an elastic junction. Assuming 
small deflections and negligible structural coupling, the three deflec- 
tion modes (flap bending, in-plane bending and torsional) are treated 
separately. 

For the bending modes, deflections occur in a plane and each rigid 
element Is constructed with three masses, m^, m 2 , m„, chosen in 
order to reproduce the same total mass, center of gravity and inertia 
of the actual blade portion this element represents. The elasticity 
of the junction is modeled by a stiffness coefficient k relating the 
local torque to the ai^le between two adjacent elements 

The equivalent blade is obtained by linking these elements and adding 
the masses at the junction 

This system is defined by the ( N + 1) coordinates of the element tips. 
The coordinate defines the position of the blade root and is either 
taken equal to zero for fixed hub conditions, or has to be defined by the 
hub equations. 

For computation purpose, it is convenient to consider that the equiva- 
lent blade is constructed with two-mass elements ( ^ central mass 

and m i ,tip mass) as sketched in fig. 4. . 


GD 005 GEN. 
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• 4. 2. 2 Definition of the blade elements 


Let C ^ be the blade linear density ( f r - ). An element of 

length £ , representing a portion of the blade between R ^ and f? ^ 

( 2 a R, ). is such that 

^ rvi, + rn^ + m^s/ 2 par-r, <totsl mass equation) 

I 

J ^Itn^ + 2^ r ir = (center of mass equation) (1) 

/ _L S yT^* Cp r *dr -rig (inertia equation) 

m2 3/ 

The values of , Tg and 1 3 are obtained from the blade mass 
tributioiuUslng the solution of system (1) gives : 


Idis- 


"^1 


I. -si^ / a. 

4 = 



c 

4 / C 

- 4 I3 / 

( 2 ) 


= 

- 1* / 4 

* 2 I 3 /£" 


Central and 

tip masses of the i th element are finally given by : 





(3) 


= 

’^■2 C^) 




4.2.3 Computation of the stiffness coefficient 

The stiffness coefficient to be associated with the hii^e between 
element ( 1 - 1 ) and L is evaluated by equatii^ the deflection energy 
for the real and equivalent blade. 

For the real blade, if and are the radial distances of 

the centers of elements ( L-1 ) and c respectively, this energy is 


W ^ 4 . 


f 


M 

El 


dr 


(4) 


p ^l O T QGA y- tvON 


G D oos GEN. 




GIRAVIONS 


Doc'--DH 2011-0. -El 

DORAND 


Page 


where M Is the bending moment, T the section inertia and E , 
the Young modulus. 


For the equivalent blade : 




= 4 


(5) 


where 0;_ is the angle between elements ( L ) and . It is 
assumed that the bending moment M is constant and equal to 6 ^. 

Combining (4) and (5) then leads to : 


- r « 


( 6 ) 


t -1 


Knowing the values of rn (equation 3) and k i (equation 6 ) 

completely defines the mathematical model equivalent to the real blade. 

The calculation of integrals I, , ^ 2 . ^3 of equation 6 

is performed simply by the trapezoidal method, Imowing the mass per 
meter and the local section inertia at a certain niunber of points. 

A sufficiently large number of points will be taken to ensure good accu- 
racy for the blade data. For example, 19 data points are required to adequate 
ly define the^section inertia and mass curves in the case of the 
DH 2011 blade (see figures 2 and 3). 

4,2.4 Dynamic equation of the equivalent blade 

The motion equation of the blade subject to forces and rotating 

at constant speed Jl is obtained by means cf Lagrai^ equations. 

See figure 4 for the notations and sign convention used. 


The potential energy U is given by 
N 


u = 4 1. k. e; 


csl 


where 0* „ V. y 

L - Oc” ® L-1 


The Idnetic energy T" is given by : 
N _ ^ N 

* 2 


T= 1 z -4 L + i r -1 ( . 


n. 

>_ 

Lsi 


1 2 2 '-' 2 ^ 2 ., 

+ T ^ . JL ) 
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Since the angles are small : 

V.- V-. _ 

Sirv 5^ = 0,,— 

enabling to written in the form : 

0. _ li _ X. (JL. ^ -±— \ (1) 

^c-1 

It is also possible to write ; 

ri- il 1^') (2, 


r- = ’"C - T 


'iT* - r- - ^ C"i- — ^- 
0 c. - N 2 ^ 2. > 


The partial derivatives of the Lagrange equations can be expressed as 
follows, using equations 1, 2 and 3 : 


k. A- / I , ^\i ^C+2 /^v 

oy,- — - -I®-’ 

--P^ = ( m. +m- ^ r'^^l 

£ 2 




, / ' 
( rn '. 4- rr> • rt 


mj rj + r;- 




't-'t 5 '*■**' ^ ^ 


The writting of these expressions may be simplified by letting : 


• f 

m./ f/ 

•^l^c -V^ 
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L 1 * ' 


:• = 21 m* r • + rn 
J=Ui 


+ mj Tj 


/f I I 

m mi+ roi.^ 


The final equation may therefore be written in the form 


e . _ k . A . u • 

ri=— ^ ^L + 1 


^c+2 ^i.4-2 ^ nr\ ♦ ^ i. rm ! X, 

0, 3 »- ^ 4 L Ot,^ 

c+t 


■**4 "'t '*• ^L<=»‘l ^Tc+1 


where 


'<^C = <-hi-'bi--,) / 

6; = .2i _ X. , U; + iilS. 

<■ 6---1 '- tc-1 

Equation 7 may be expressed in the following matrix form 


EM,Z-«-£K.Z. -I- Jl^ EC. X (8) 


where 


Fi 





♦ 

^ a 


«* 

z = 

• • 



5n 


• • 
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F being the column vector of the forces and that of the displace- 
ments. 


EM , £K. and £c are square matrices ( N x M ) which represent 
successively the effects of mass^ elasticity and centrifugal forces. 

The elements of these matrices may be calculated from the following 
formulas { M ^ being the term in the i th line and j th column). 


Matrix 

EM 


£M. 


-i- m; 

4 ^ 


E = no- -+ ± r rrli 


Matrix* 

EK 


E 

s liL + k- af 4- 


a .2. , 


— - k (^*1 -v k L-va L - 


EKr = 


& ' p ^ 

^ L-hT. 


Matrix 


_ ^i±£L 


E<^.- 


EC.'-*' , 






Remark : The non-defined terms of these matrices are zero. 

It may also be seen that these three matrices are symmetrical (since 

+ ^L + 1 ^ which simplifies calculation. 

Equations 9, 10 and 11 thus enable the blade motion equation to be ex- 
pressed in linear differential matrix form, which relates the flapping 
displacements of the N points selected on the blade (at the junctions 
of the equivalent blade elements) to the aerodynamic forces applied to 
these points. 


F=GM.X <4- [EK+JI^ecJ.Z 
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This equation is easily resolved with the Fourier coefficients of the 
displacement forces, since this amounts to a matrix inversion. 

Thus we can write : 

+ £K -F GC ] ^n. 

where is a vector containing the n-order Fourier components 
(sine or cosine) of the forces on the various points, 

same definition for the displacements 
The n-th harmonic is thus given by : 

[_ -F EK + ec]"^ 

Remark : This equation should be applied twice for each harmonic, 
once for the sine and once for the cosine. 

4.2.5 Calculation of Flapping natural Vibration Modes 

Takii^ the basic matrix equation : I 

r « E H .Z + [e« f EC ] Z 

it is seen that the resonant frequencies can be easily calculated by 
writing that for these frequencies this equation cannot be inverted. 

Let to„ be a resonant pulse , giving : 

[-cu^ en -h EK + ] Z 

The conditions cf resonance impose : 

determinant EM, EK J «0 

z. 

This equation is true if co is a eigenvalue of 
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The various natural frequencies are obtained by calculating the eigen- 
values of EM-"* QEK + JI^EcJ. 

It is apparent that the number of values obtained is finite and equal 
to the number N of elements of the equivalent blade, whilst in 
reality it is infinite. This raises a problem of accuracy and therefore 
of convergence of the method. 

We have, for example, evaluated the rate of convergence by calcula- 
ting the natural frequencies of the first three modes of vibration of 
a uniform bar, for different numbers of elements constituting the 
equivalent bar. Two cases of calculation were considered : zero speed 
of rotation and a rotation of 33.3 rd /s. 

The results are shown in figures 5 and 6. 

Since we theoretically know the required values for a uniform bar, 
we note that : 

- the method provides good accuracy for the first 3 modes with a small 
number (approximately 6) of elements- 

- convergence is more rapid for zero rotational speed than for 

= 33. 3 rd /s. 

The application of this method to the DH 2011 blade enabled us to 
determine the first three modes (flapping and the first 2 bending modes) 
and to compare them with the values calculated previously. The calcu- 
lation was made for a 6-element equivalent blade. See figure 7. 

4.2.6 Calculation of blade bending stresses 

The blade motion equation enables us to know the Fourier coefficients 
of blade distortion. 

Knowing the cross-section thickness of the blade spar is therefore 
sufficient for determing local stress in this spar. 

See figure 8 for the followii^ calculation. 

If is the ha If -width of the spar, we can express the stress 

as : 
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For small angles, however, ^ « Vv,6^ 


whence 

^1 = 

£6- . 

^ ti 

Since 

- 

V- - ^ 


then finally : 



This equation can be written in matrix form : 

G" es.z 

where CT is a column vector containing the stresses at the M points 
on the blade. 

This equation can also be written using Fourier coefficients : 


cr^ = e s , 

giving finally the Fourier coefficients of the stresses. 
Matrix £S is defined as follows : 



6.K- 

1 

ei eui 

L-1 

ES. e 

_ e,k. 

* ec 

ES"" = 

e .hi 

1 




The stresses can therefore be easily calculated by simple matrix pro- 
ducts. 

4. 2. 7 Calculation of hub vibrations 

Vertical hub vibrations are determined by calculating the blade hinge 
force. In order to obtain the total force, both blades must be considered, 
which results in doubling the even- order Fourier coefficients and eli- 
minating the odd-order coefficients. Odd-order vibration harmonics 


G D 005 G£N. 




GIRAVIONS 

DORAND 


Doc'. _DH.20.11-C^E1 

Page “18- — 

cannot occur in a two-bladed rotor. 

The force applied by one blade is written : 

P = e, (1) 

where is the centrifugal force. 

Equation 1 can also be written in the form : 


F ^ 



e, 




Since the Fourier coefficients of and Ja are known; this 
equation gives the coefficients of the force applied to the hub by the 
blade. 

The constant term corresponds to lift. 


5. DESCRIPTION AND USE OF THE FORTRAN IV PROGRAMS 

As already stated, the computational program consists in fact of three indepen- 
dent programs : 

BFHA : Blade forces harmonic analysis 

BSTM : Blade stresses tranfer matrices 

STRESS : Computation of blade deflections and stresses. 

The following is a brief description of each of these programs, for which the 
listings are given in Appendix 1. 

5. 1 BFHA ; Blade forces harmonic analysis 

The distribution of forces along the blade is given at a certain number of 
points on the blade (17 in the present case) and for different azimuthal 
positions (24 in the present case, i. e. every 15”). The purpose of the pro- 
gram is to provide the Fourier coefficients of these forces (at the same 
radial points as the data). 

PROGRAM INPUTS : Read operations are performed in the following order : 

commentary cards : the number of cards is not limited (format 20 A 4). 

- ^ INPUT DATA : this card (punched starting with column 1) indicates that 
commentary is complete and that computing data will follow. 
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the next card contains in order : 

. simulation identification n* in F 6. 2 format (column 1), 

number of azimuthal positions where the forces are given, in 
13 format (column 26), 

. number of radial points in 13 format (column 37), 

, number of harmonics required in 13 format (column 47), 

unit (Newtons per meter in present case). This enables a factor to 
be given for changing from the forces as given by the simulations 
to forces in Newtons/meter. NW/M is contained in column 68. 

the next cards define the positions of the radial points (values in frac- 
tions of radius) measured from the blade root. Format F 10.4 is 
used. Seven values are given per card (columns 11, 21, 31,41, 51, 

61 and 71). 

Set of cards corresponding to one simulation (from EMC program). Put 
the last 4 cards of this set of cards ahead. 

Card ^ END (column 1) indicates there are no more cases to be 
treated. 

If this card is not present, the program starts reading new data (same 
presentation as before : commentary cards, etc. ). 

PROGRAM OUTPUTS : 

line-printer output containlt^ input data, and the Fourier coefficients 
harmonic by harmonic for all the radial points defined by the input data. 

punched-card output containing the same data and directly usable for 
computing stresses. 

This set of cards is entitled PUNCH 1. 

5. 2 BSTM ; Blade stresses transfer matrices 

The purpose of this program is to compute the transfer matrices EM, EC, 
EK and ES defined in § 4. 1. 2 and the natural vibration modes. 

Since there are a greater number of Inputs and outputs for this program, 
the following description refers to A^-penHix 1, which gives the listii^s. 

PROGRAM INPUTS 

1) - Two cards containing : blade identification, rotor radius, distance 

between the axis of rotation and the blade root, unit of lei^h, 


PKOTO05A V - (Yprs 
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material modulus of elasticity and its unit, and the number of 
data points of the blade characteristics (FORMAT 501). 

the next cards contain the radial positions of the blade characte- 
ristic data points (FORMAT 502). 

the next cards contain the blade longeron thickness (FORMAT 502). 
card containli^ the units of mass and Inertia (FORMAT 500). 
option 1 card : 

Two possibilities : COMPUTED BLADE ELEMENT PARAMETERS 
or GIVEN BLADE ELEMENT PARAMETERS 

option 2 card : 

Three possibilities : CONSTANT VALUE OF ELEMENT LENGTHS 

GIVEN VALUES OF ELEMENT LENGTHS 
OPTIMIZED VALUES OF ELEMENT LENGTHS 

if option 1 is "COMPUTED", the next cards contain the mass per 
meter (density) and the cards contain the blade section inertia 
(FORMAT 502). 

card NBE containing the number of equivalent blade elements 
(FORMAT 503) 

if option 2 is "GIVEN", the blade element lengths are read 
(FORMAT 502) 

if option 1 is "GIVEN", the spring constant^ and the masses 
assumed to be concentrated at the center and extremity of each 
element are read (FORMAT 504) 

option 3 

Two possibilities : PRINT BLADE ELEMENT PARAMETERS 
NO PRINT 

option 4 

Three possibilities : PUNCH AND PRINT TRANSFER MATRICES 

PRINT ONLY 
NO PRINT NO PUNCH 

option 5 

Two possibilities : CONTINUE 

EIGENMODES 

If option 5 is "EIGENMODES", the next card gives the speed of 
rotation in radians/second (FORMAT 505), and option 5 is reread. 
This loop repeats until a "CONTINUE" card is encountered. 

option 6 

Three possibilities : NEW BLADE CONFIGURATION 

CHANGE VALUE OF NBE 
END OF COMPUTATION 
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= if option 6 is "NEW BLADE CONFIGURATION”, the sequence 
returns to 1. 

= if option 6 is "CHANGE VALUE OF NBE", the sequence returns 
to 2. 

NOTE : Cards or groups of cards marked in the margin by are 

necessary, whilst those marked by"=" depend on the contents 
of certa in options. 

PROGRAM OUTPUTS 

The line-printer output is shown in figures 9 through 12. 

the first page lists the mechanical characteristics of the blade :spar tick- 
ness, mass distribution and section inertia, followed by the total mass 
average section inertia and the 1st resonant frequency of a uniform 
blade having the same average characteristics. 

the second page lists the characteristics of the equivalent blade, i. e. 
the masses constituting each element and the spring constants Unking 
them. The four matrices EC, EK, EM and ES are then printed out. 

the third and fourth pages gives two examples of natural mode compu- 
tations (EIGENMODE option) for rotational speeds of zero and 318 rpm. 

The punched-card output contains matrices EC, EK, EM and ES, in addition 
to certain characteristics such as radius, the number ctf elements, etc. 

This set of cards may be directly used for the third program and is 
entitled PUNCH 2. 

It should be noted a punched-card output is obtained only if the PUNCH 
option has been requested. 

5. 3 STRESS : Stress and vibration computation 

The purpose of this program is to compute the rotor stresses and vibration 
due to the applied forces (output of program BFHA), using the matrices 
defining the equivalent blade (output of program BSTM). 

PROGRAM INPUTS 

The input cards are stacked in the following order : 

PUNCH 2 (cards produced by BSTM) 
a punch option card : 
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PUNCH in column 1 


PUNCH in column 40 


unpunched card 


punched-card output of the blade 
distortion Fourier coefficients 

punched-card output of the blade 
distortion as a function of blade 
azimuth ( if'' ) 

neither option selected 


PUNCH 1 (cards produced by BFHA) 

Three cases are possible following this computation : 

1) There are no more computations and computation therefore stops : 
add the following itput cards : 

1 unpunched card 

1 mt END carr’ (column 1) 

2 unpunched cards 

2) It is required to enter another computation which uses the same matrix 
data for the equivalent blade, only the system of applied forces being 
changed. 

In this case add the following ii 4 >ut cards : 

1 punching option card 

the new stack of PUNCH 1 cards produced by BFHA followed by the 
same procedure as before. 

3) It is required to enter another computation with changing : the matrix 
data for the equivalent blade and the system of applied forces. 

In this case add following input card 

1 unpunched card 
1 ^ END card 

the new stack of PUNCH 2 cards produced by BSTM. 

Followed by the same procedure as before. 

PROGRAM OUTPUTS 

The following data are printed out in succession ; 


EC, EK, EM and ES matrices 

Fourier coefficients of the forces and coefficients C,j. and M.j. (force 
and moment) for the rotor 

Fourier coefficients of the blade distortion 
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blade distortion as a function of . 

Fourier coefficients of blade stresses 
^ stresses as a function of 

^ Fourier coefficients of vibration (vertical force) 

Fourier coefficients cf stresses at 0. 45 K 
Fourier coefficients of stresses at 0. 7 R 
^ vibration as a function of ^ 

^ stresses at 0.45 R as a function of . 

various parameters of rotor operation with its cyclic control in parti- 
cular (including multicyclic) 

Results marked ^ are punched out onto carts independently of the PUNCH 
option. 

Remark 

Certain results obtained on cards, such as the 0.45 R and 0. 7 R stress 
Fourier coefficients and vibration, have the same format as that used 
for the input data of the multicyclic analysis program (see reference 4). 

5.4 Conclusions concerning presentation of results 

This report does not cover all results obtained on the line printer, but 
tables have been prepared comparii^ real cases and simulations for the 
stress Fourier coefficients. 

In addition curves, obtained by means of a curve plottii^ program and 
comparing stress curves for real and simulated cases, are presented. 

SUMMARY OF RESULTS 


The following presents the results of the stress calculations at 0. 45R for the 
15 simulations at ^u = 0.4 performed with the Evans -Me Cloud program. 

These 15 cases were computed with an equivalent blade consisting of 6 rigid 
elements of equal lei^ht. It has been estimated from the natural frequency calcu- 
lation that 6 elements provide sufficient accuracy for computing blade stresses. 

The Fourier coefficients of the computed stresses are given in Table I and sh 
be compared with the measured stresses given in Table IIo; In order to facilitate 
comparison, measured and computed stresses have been plotted against ^ on the 
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same curve for each of the 15 cases. 

Examination of the results leads to a first remark concerning the radial dis- 
tribution of the blade stresses. The stresses measured durii^ the runs are 
maximum at 0. 45R. For example, for the 15 runs analyzed, the mean stress 
is approximately 7.75 hectobars at 0. 45R andS.Ohectobars at 0. 7R, The computed 
stresses, however, show a maximum at 0. 7R. For example, the mean stress <3o 
is approximately 6L0 hectobars at 0. 7R and 4.8 hectobars at 0.45R. (See Figure 13, 
where these distributions have been plotted for Run-PT 9. 03). 

Since total blade lift is the same in either case (measured and simulated), 
there are two main reasons for this anomaly : 

- the distributions of aerodynamic loading are certainly different for the real 
and simulated cases, 

- the mechanical characteristics of the blade, i. e. section density and inertia 
are not accurately known. 

This remark led to plotting the measured and computed stresses as ratios of 
their mean values in order to be able to compare the two curves. 

Comparison of computed and measured results does not allow positive conclusioJi 
to be drawn. It is possible, however, to state that the simulated cases contain ver 5 ?[ 
small 3rd and 4th harmonic components, whilst on the contrary the real cases 
contain high values. 

Correlation^for the 1st and 2nd harmonics is much better, as shown by 
Tables I j and Hi, in which the relative amplitudes and signs of the coefficients 
are very similar. For example, it is interesting to note the variations of the 1st 
and 2nd harmonic coefficients for runs 14-10 to 14-13, i. e. decreased 1st harmonic] 
and increased 2nd harmonic. The overall correlation remains nevertheless rather 
poor, as shown in Figures 14 to 28, because of the absence of higher harmonics. 

When this was noted, simulation with a smaller airfoil stalling angle was 
requested (12'* instead of 16°) in order to attempt to introduce larger lift discon- 
tinuities, which would therefore contain greater values of higher order harmonics 
of aerodynamic force. This charge did not produce the expected result, the 3rd 
and 4th harmonic contents remainit^ far too small. This is explained that the 
stalled sector for the blade at ja. = 0.4 occupies only a very small part of the 
disk. The value of stalling angle was therefore returned to 16“ . 


7. CONCLUSION 

The puJTJOse of the work described in this report was to create a program for 
computing helicopter blade stresses from a knowledge of the local forces applied 
to the rotor disk. 
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Local forces can be computed by means of the Evans - Me Cloud simulation 
program for the DH 2011 rotor. 

During the DH 2011 rotor runs in the 40' x 80’ Ames windtunnel in March 1971 
the blade bending stresses were measured for various aerodynamic configurations of 
the rotor (in particular, multicyclic control). Fifteen of these configurations were 
selected and simulated with the Evans - Me Cloud program. 

CD's stress program enabled blade bending stresses to be computed for 
these 15 simulations, the purpose beir^ to examine the correlation between mea- 
sured and computed values. 

Comparison showed fair agreement for the 1st and 2nd harmonics of stress, 
but very poor correlation for higher order harmonics (the computed values 
containing only very small amplitudes). 

Consequently, although the Evans - Me Cloud program produces good results 
for rotor performance, it cannot be used for analyzing dynamic phenomena, such 
as vibration or stresses. 

In order to obtain a better description of blade dynamic loads, several 
effects need to be included in the simulation program : 

- unsteady aerodynamics 

- blade flexibility 

- wake effects 

- Mach number effects. 

It is, however, difficult to state which of these effects has the greatest influen- 
ce on the 3P and 4P harmonic contents of local aerodynamic loads. 

Notwithstanding the inclusion of these effects in the simulation program would 
certainly improve the accuracy erf aerodynamic force computation. 
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This appendix contains the listir^s of the FORTRAN IV programs used for the 
stress computation- 

The three main programs are : 

1) BFHA : Blade forces harmonic analysis ; which requires the subroutine 
FHARM3. 

2) BSTM : Blade stresses transfer matrices ; which requires the subroutiaei 
ATEIG, B2C, CLEAR, EIGMOD, HSBG, IPOLl, LICOM, MINV, MOMEN' 
PRODl and STRAM. 

3) STRESS : Stress computation. 

This program is composed in fact, of a short MAIN program and of a 
subroutine named STRESS. The other subroutines required are B2C, 
FSUM3, IPOLl, LICOM, MINV, MOMENT and PROD. 

The listings of the subroutines are also joined to this appendix. 
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BLADE FORCES HARMONIC ANALYSIS 

DIMENSION H<2) «TYPE(2» »UF(2)»T<20I fRI20»*OeLTAI9ltCOMI7> 

COMMON F0<20I.FN< 12«20t2) ,FI24*20f 21 

DATA H.TYPEtDATIN,END/»A%*BS*FZ*,*FX*,**INP»,*»ENO»/ 

500 FORMAT(20A<tI 

501 FORMAT ( F6,2«9X. *RPM= S F6. 1» *NPTS=* 1 1 3* 2X, • NPOST=* »I3f 1 X» 

1*NHARM=' ,I3aX,*FACT0R=»* ElO*3*2A4/( • X/R = '*7F10,4n 

502 FORMAT! I3f6ElU4I 

503 FORMAT! 7 A4, 1 2X, F5 • 2t 19 X, F5, 2/lOX# 9F ?• 2// 

1 3X*F6,2*6X,F7.2*8X.F6.3t6X,Fl0.7) 

600 FORMAT! »1SET NUMBER' ♦13/1 

601 F0RMAT!10X,20A4» 

602 FORMAT!//' RUN I DENTIFICATI ON* , 2X . F6.2 * 5X» ' RPM= ' ♦ F6*l// 

IIX, 13. 'POINTS PER STATION' . no, • STATIONS AT X/P =' //! I X, 18F7.4H 

603 FORMAT! 1X,A2, I3,17F7.41 

612 FORMAT! /• FORCES VALUES MUST BE MULTIPLIED BY',1PEU,3, 

1* TO OBTAIN*, IX, 2A4//' AMES DATA CAROS INPUT •••*/) 

604 FORMAT! * IFOURIER COEFFICIENTS FOR ',A2,I10,' HARMONICS ANALYSED*/! 

605 FORMAT! 1X,A2,1X,A1,I2,4X, IP 1 OEl 1. 3 /! 1 1 X , lOEl 1, 3 J ) 

700 F0RMAT(60X, 'OPTION CARD 

701 FORMAT! F6.2 ,4X,* AL PHAS = * ,F6 ,2 , ' C JR= • , FI 0,7 , ' MU=',F6,3,' SL=», 

1 F5.2,' THETA=' ,F5. 2/F6, 2,'OELT' ,9F7, 31 

705 F0RMAT!A2,1X, II, I3.3X, lP7E10.3/a0X,7E10.3n 
NSET=l 

6600 WR ITF! 6,600 INSET 
WRIT6I7,700I 
5500 READ! 5,500)T 
WRITE(6,601 IT 
y<RITE!7,500)T 
IF(T! 1 J.EO.ENOI • STOP 
IF(T! n.NE,DATIN 1 GO TO 5500 

READ! 5, 501) RUN, R PM, NPTS, NPCST , NHARM, F ACTOR, UF, CRd), I=l,NPOST) 

READ! 5, 503) COM, SL, THETA, DELTA, ALF AS, OMR, AMU, C JR 

WRITE!6,601) COM 

RAD = 19.68 

OMEGA = OMR/RAO 

RPM = 9.549*0MEGA 

FACTOR = 0.962540HeGA4424!RAD40.3048»**3 

WRITE! 7,501 ) RUN , RPM, NPTS, NPOST, NHARM, FACTOR, UF, I R! I ) ,I=1,NP0ST) 
WRITE!6,602) RUN , RPM, NPTS , NPOST ,! R !I », 1= I , NPOST ) 

WRITE!6,612) FACTOR, UF 
NC=INP0ST-1 1/6+1 

1001 DOl K=1 ,NPTS 
11=1 

12=6 

1002 002 J=1,NC 

1003 D03 NTYPE=1,2 ' 

3 READ!5,502) KPS I , IF ! K, I , NTYPEI , 1= 1 1, 1 21 

n = I 2 +l 
12=12+6 

IF! I2.GT,NP0STI I2=NP0ST 
2 CONTI NUE 
1 CONTINUE 

WRITE 1 6, 603 ) ! (TYPE! NTYPE I , Kf!F(K, 1 ,NTYPE) , 1=1, NPOST I ,K=1,NPTS) , 
1 NTYPE = 1,2) 

C START THE HARMONIC ANALYSIS 

NTYPE = I 

C STATEMENT 1004 SUPPRESSED FOR ANALYSING FZ ONLY CNTYPE * 1) 

Cl 004 00 4 NTYPE-1,2 
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L=1 

N=0 • 

CALL FHAPM3 (NPTStNHARf^, NPGST.NTyPE ) 

WRlTE{6«604i TYPE < NTYP6 > , NHARM 

WRITE <6, 6051 TYPE (NTYPE } , H( L J *N* « F0( K J , K= 1, NPOST) 
WRITE<7»705> TYPE t NTYPE J »L tN t< f 0 ( K > ♦ K=l ♦ NPOST I 
2004 00 4 NaltNHARM 
3004 DO 4 L=1.2 

WRI TE( 7*705) TYPE(NTYPE) *L* N* ( FNt N ,K * L » » K=1 , NPCST I 
4 WRITEI6,605I TYPE < N TYPE ), H< LI , N* < FNC N* K * L )» K=1 * NPOST) 
WRITE(7,70l) RUN. ALFAS*CJR*APU,SL*THETA, RUN, DELTA 
NSET*NSET+1 
GO TO 6600 
END 
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C BLAOe STRESSES TRAVJFER H4TRICES PROGRAM S 

C0MM3M/BL A0E/BEL! 20 K 3ERC 20 > t SECM* 201 , BETMJ 20) tBEKI 2 01 t 
1 A(2O),8<20) 

COMMON/TRAMS/OMEGAf E SI SO) tEM160i fECtSO) tEKClOO I 
C0RM3N/«0R</Ml( 20 Iti<2l20)» W3I20) *«4I 400) t W5I 4001 
DIMENSION TEXTI19) «C0DEI7)f BDIBOIf BII30)vHEI30)»RS(30l» 

I UU3),UM<3)*UH3)*Uei3)»BUtOI4) 

DATA C0DE/*GIVE* f •0PTI*t*PRlN*f "EIGE't'NEH • t * PUNC • t *CHAN» / 

500 FORMA Tl 2044) 

501 FORM ATI • IBLADE » 1 4A 4, 7K, • R= • tlPEl 0. 3* IX • »R0«* • El 0. 3 #l X*3 A4 /I Xt i 

l» £*• *E10,3* IX*3A4*4X vUMBER OF BL40E DATA POl NTS't 13 ) 1 

502 FORMAT! 10X*IP7EI 0.3) I 

503 FORMAT! 'IMJMBER OF BLADE ELEMENTS «S13) 

504 FORMA TH0X,3E10.3) : 

505 FORMAT! ‘OMEGAaSPS. 1) ' 

506 F0RMATIA4»lXt lU : 

600 FORMAT! /• *PROGRAM OPTION* *.2044) j 

601 FORMAT!/' R AO IDS IPE 12, 3. I X,3A4/' TOTAL MA SS* .El 2 . 3 *l X. 3 44 / 

I' AVERAGED SECTIONAL I NERT lA*. E15.3. IX. 3A4/ 

2* UNIFORM BLADE FREQUENCY* ,E12. 3 ♦' RAD/SEC*) 

S02 FORMAT! » EL EMENT* . 7X. • LENGTH*. 3X. «R AD. DI ST, SPRI NO CTE. CENTRAL 
IMASS TIP MASS*//(4X.I3.4X.IP5E12.3)) 

603 FORMAT!/* BLADE STRUCTURAL DATA*/* EL EMENT* ,7X . * RAO. 01 ST* ,4X . * THI C 
I KNESS*,2X. * MASS KG/ M)*,2X.* SECT. INERT *//!4X, 13. 4X.1P4E 12. 3) ) 

C 

c 

1 N3PT « I 


READING OF THE FOLLOWING OATAS= 

BLID=BLADE IDENTIFICATION 

BRA0=8LA0E RADI US. RO*OIST ANCE FROM THE ROTOR AXIS TO THE/ 
BEGINNING OF THE BLAOE, 

ULaLENGTH UNITS. 6«ELAST IC ITY MODULE. UE=MODUL E UNITS. 
NR»NUMBER OF BLAOE DATA POINTS. 

RE AO (5. 501) ! SLID! I ) ,1*1 , 4 > . BRA D. RO. JL . E , J E. NR 
WRlTE!6.50l) IBLIO! n .1*1.4) .BRAO.RO .JL.E.UE. NR 
E=U. E+07)*E 

READING OF RS AND HE. 

HE*THIC<NESS of the BLAOE LONGERON. 

RS *RAOIAL DISTANCES WHERE BLADE CHARACTERISTICS 
ARE MEASURED, . 

REA015.502) ! R$! I) « I *l • NR) 

REA0I 5.502I ! HE I I) . I* l.NR) 

READING OF HASS UNITS UMI AND INERTIA UNITS !UI). 

REA015.500) UM.Ut 

READING OF OPTION l« 

TWO CASES COMPUTED BLADE ELEMENTS PARAMETERS 

GIVEN BLAOE ELEMENTS PARAMETERS 

REA0(5.500) dPTl.TEXT 
WRITE 16.600) OPTl.TEXT 

•■■ ■"( 1 

READING OF OPTION 2= 

THREE CASES CONSTANT VALUE OF ELEMENT LENGTHS 
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VALUE OF B.ENENT LE*4GTHS 


OPTlHIZeO VALUE OF ELEMENT LENGTHS 


E 1 


REAOC5.500I 0PT2,TEXT 
WRITE I 6, 6001 0PT2*TEXT 


RESET TO ZERO IS.P* CLEAR > 


CALL a EAftI 1 » 

IF OPTION I NOT EQUAL ‘GIVEN* READ BLADE DENSITY ( BDI AND 
SECTION INERTIA IBII. 


IFIOPTl.EO.COOEUH GO TO 2 
ReA0(5.502l ( 30 m«I«l,NR) 

REAO(5,502) I BI (I ) • I =l t NR I 

WRITE! 6f 603 » I I . R S( II ♦ HE ( 1 1 # BOI U t BU II » I « l» NR > 
N8E=E0UIVALENT BLADE NUMBER OF ELEMENTS, 


2 REACH 5, 5031 N8E 

RESET TO ZERO. 

CALL CLEAR! 21 

[F(OPT2,EO.COOE!l II GO TO 5502 

IF OPTION 2 EOUAL»GI VEN'READ THE LENGTHS OF BLADE ELEMENTS, 


IF!OPT2*E3.CODE! 2 » I S3 TO 30 

IF OPTION 2 NOT EOJAL'GIVEN* OR • OPT I • , COMP JT E THE LENGTH 
OF THE BLADE ELEMENTS, 

BELS * !BRA0-R01 ZFL0ATIN8EI 
1003 00 3 1=1, NBE 

3 BEL! I I s BELS 
GO TO A 

5502 READ! 5,502I ! BEL 1 1 1 , I « 1, NBE I 

IF OPTION 1 NOT EQUAL 'GIVEN* COMPUTE THE MASS, GRAVITY 
CENTER AND INERTIA OF EACH BLADE ELEMENTS, 

4 IFIOPTl,EO,COOE!l)I GO TO 20 

COMPUTE INVERSE VALUES OF SECTION INERTIA AND STORE IN W2 
ID05 00 5 1=1, NR 

5 H2(II=1./B1(I1 

COMPUTE BLADE ELEMENT TOTAL HASS, CG AND INERTIA |SM,SMl,SM2I 
ASSUME FREE HINGE BLADE !8EK!II=3,I 

BEK! 1 1=0, 

XI = RO 
Kl = I 

XMl=RO+BEL!l I /2, 

KMl=l - 




c L80P FRO« 1306 T3 6 COMPUTES 8LA0E ELEMENTS PARAMETERS® 

C BECM*CENTRAL MASS 

C BETM=TIP MASS 

C BE< ^SPRING CONSTANT. 

C 

C 

1006 DO 6 1=1. NBE 

8ERU > a XI ♦ BEUn 
X2 « BER«n 

CALL MOMENT(Mlt2.RS*XI.X2« Xl.80tKl.NR.K0) 

SM = wmi 

SMI = Wl( 21 /BELf I) 

SM2 = WU3)/BEL{n*>^2 
C 

c ■ 

c COMPUTE THE EOJIVW.ENT MASSES I EMI. EM2» EM3I 

C 

EMI = SH - a.’i'SMl ♦ 2.*SM2 
EM2 = 4.#<SHl-SM2) 

EM3 = 2.*SM2 - SMI 
C 
C 

C COMPUTE BLAOE ELEMENTS CENTRAL AND TIP HASSES( BECM.BETM) 

C 

BECMtl) = EM2 
0ETMI I ) = EM3 
IF4 1.60.1) GO TO 7 

BETMU-ll = 6£TM(I-ll ♦ EMI 
C 

C COMPUTE BLAOE ELEMENT SPRING CONSTANT BEK 

C 

XM2=BERU-l) + BeLCn /2. 

CALL MOMENT! H 1,0, RS. XM 1 , XM2 . XMl , H2.KM 1, NR.KO) 

BEKn) = 6/NlIl) 

XM1=XM2 
GO TO 6 
7 BETMO = 6K3 
6 XI = X2 
GO TO 40 
C 
C 

C IF OPTION EQUAL 'GIVEN*. 

C READ GIVEN VALUES OF BLADE ELEMENTS PARAMETERS 

C 

20 BRAOaRO 

1021 00 21 J=1,NBE 

REAOI 5,504) BE KI J| .BECMI J ) ,BETM( J) 

BER(J) = BRAD ♦ 8EU Jl 

21 BRAD = BERIJ) 

C 

c 

c COMPUTE AVERAGED PARAMETERS 

C 

GO TO 40 
30 CONTINUE 

C . ■ 

C CALL VARBELINR.NBE) IN CASE OF QPTIMUATION OF THE LENGTHS 

C ^ . 

40 IFINOPT.EO.O) GO TO 41 

BRASS * 0. . 
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BftAQ » BEMNBE) 

8 X A V “ 0* 

1050 00 50 1=1, NBE 

8M4SS = 8M4SS ♦ 8ECM(II ♦ 9ETMC II 
IFf I.EO.l) SO T3 50 

BI AV=BIAV+BEKm*CBELU-l l+BEU II »/2, 
50 CONTINUE 

BIAV»8IAV/{ FL0ATINBE-1I*EI 
BMASS * BHASS ♦ BEXNO 


COMPUTE THE EOJIVALENT UNIFORM BLADE El 5ENFREQUENC T 

□ MU6 = IS.AfrSORTIEfBI AV/IBMASS««<8RAD-R0l**3n 
WRITEI6,6011 BRAOfULf BMASS, UM, 8 1AV,UI , OMUB 


READ THE OPTION 3. 

TWO POSSIBILITIES PRINT BLADE ELEMENTSS PARAMETERS* 

NO PRINT. 

41 WRITEI6,503) NBE 
RFA0(5,500> OPT, TEXT 
IF(0PT.NE.C00£C3 n G3 TO 200 

WRITE(6,630) OPT, TEXT 

WRITE 16,602) ( I, BEU I) ,3ER( I ),3EK< I) ,BECM< I) ,BETMI I ) ,1=1 , NBE) 
READING OF THE OPTION 4 

THREE CASES PUNCH THE TRANSFER .MATRICES 

PRINT ONLT 
NO PRINT OR PJNCH 

200 READ(5,500I OPT, TEXT 
WRITEI6,600) OPT, TEXT 
IS = 0 

IFIQPT.EO.CODEU) ) IS=6 
IFItlPT.EO.COOe(6) I IS=7 


COMPUTATION OF THE BLADE STRESSES TRANSFER MATRICES EM, EX, EC 
CALL STRAM(NBE,NR,IS,RO,BLID, HE,RS,EM,ES,EC, EK ) 


READ THE OPTION 5 CONTINUE 

EI3ENMOOES 

207 REA0I5,500) OPT, TEXT 
MRlTEt6,600) OPT, TEXT 

IF OPTION 5 EQUAL 'EIGEN* , READ THE VALUE OF OMEGA , 
COMPUTE EIGEN MODES OF THE BLAOE,ANO READ OPTION 5 AGAIN. :; 

IFIOPT.NE.CODE(A)) GO TO 210 

REA0<5,505I OMEGA 

CALL EIGHDDIEM, EC, EK, OMEGA, OMUB, NBE) . 

GO TO 207 

READ OPTION 6 oT 
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210 READ< 5,500J OPT, TEXT 
«RITEI6,600I OPT, TEXT 
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THREE CASES NEW BLADE COMPUTATION 

:HANGE THE VALUE OF MBE 
END THE COMPUTATIONS 

IFOPT.EO.CODEISn GO TO 1 

NOPT = 0 

If (OPT.EO.CQOeiTn GO TO 2 

STOP 

END 
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COMMON CF.BB^RO 

DIMENSION BERI15>,ECJ 15. 3) . EKU 5. 5} .EMU5* 31 • E SI 1 5 ,3 ) «T (301 »Rf 23 ) . 
102(21 .101(21, IMl (15),I»I2(15),SIGM&( 15. 2 1 » , F 2IU I 36, 31 ^FC( 3,211 . 
2W1I20I ,W2(20) .W3I20) ,M^( 201 ,^(225 (.<(6(225 ),W7( 225), F2RI8, 15), 

3BEL ( 15). 1021 2),SLID( 41 

5700 REA0(5.700) BLIO.N BE, CF. BB. RO. ( BER ( I ) . I« 1. NBE ) 

700 FORMATCIBLAOE I OEN TIFIC A TI ON • ,4A4 ,2 X , • N8E=* , 13, 3X, • CF=* . 
IIPE10*3,« B6»*,E10.3/( 10X,7E10,3) ) 

URITE(6,700) BUO.NBE.CP.BB.RO, IBERII), I«1.NBE) | 

|F(N8E,E0,0» STOP I 

NBE2 = NBE«*2 i 

CALL STRESS(N6E,NBE2,BER.EC,EK,EM,ES. BEL, IWl. IW2.SIGMA.F2R. 

IWl. W2. W3,M4,W5,U6,W7) 

GO TO 5700 j 

• END ! 
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SUBROUT I^^E STRESSJN8E*N8E2*BER*eC.EK*EM*ES.BELfIWl ,I H2 t SI GMA,F2R, 
I Ml W9tW6«W7) 

COI4HON CFtBBtRO 

OMENS ION BERCNBE I, ECI NBEf 3 )t E<J NBE, 5) , ENI N8E* 31 1 ESt N9E t3 ) » 

I7f20) «R(20t «02I2) 1 1 OU2 > « Ml CNBEI t M2 IN BE) « S IGMAI NBEf 2 1 It 
2F2R1I36«3) tFCI3*21) tUll I I *W2I I ItHBI I )*W^I 1 )»W5(N8E2», 
3W6(NaE2)tW7(NBE2)tE2R<6tNBE)«BELtNBE It ID212) 

DATA OATINtENOt!02» YES /■ •INP* ♦END* t* A« t *8 • t • PJNC« / 

599 FORNATI* X/R= • * lOF 7. 4/ C 1 OX, I 0F7, 41 I 

600 FORMATdXt* Z * • t A1 , 12, 2 X, 10F7 *3/ ( lOX , 10F7, 31) 

601 FORMATI* PSI «• ,F5, 0, • DEGREE S« / 1 lOX ,10F7,3 1 ) 

602 FORNATI'l*) 

603 FORMAT!//* DIRECT INTEGRATION OF FORCES*/* 10004CT«* ,F 7,3t 

1 • 1000*HT**tF7, 3) 

:604 FORMAT!*! BLADE DEFLECTION I100*Z/R) RADIAL DISTRIBUTION FOR OIF 
lERENT AZIMUTHAL ANGLES*) 

605 FORMAT! lXt20A4) 

606 FORHAT!IXt2A4tIX,lP7El0.3/! 10X,7£10,3) ) 

607 FORMAT!* 1 BLADE DEFLECTION !100*Z/R) FOURIER COEFFICIENTS*) 

610 FORMAT!*! STRESS ! 10000*SI GMA /E I F3URI ER COEFFICIENTS*) 

611 FORMAT!*! STRESS !10000*S I GMA/E I RADIAL DISTRIBUTION FOR OIFFERE' 
IT AZIMUTHAL ANGLE S* ) 

620 FQRMAT!IX,*SIG * t A 1 1 1 2t 2X, lOF 7, 3/ ! lOXt lOF 7, 3)) 

701 FORMAT!* EC* t7X tl P3 ElO, 31 

702 FORMAT!* EK * t 7Xt 1P5E 10, 3) 

703 FORMAT!* EM* t 7Xt IP3 E10.3 ) 

704 FORMAT!* E S* 1 7X ,l P3 ElO. 3 ) 

705 FORMAT! 20A4) 

706 FORMAT!2A4t2X,lP7£10.3/!10X,7El0.3 ) I 

708 F0RMAT!F6.2t4Xt* RPM=* tF6. 1 1 • NPTS=* ,13 tZXt* NPOST=* , 13, IX,* NHARM=», 

1 I3,IX,« FACTORS*, E10.3,2A4/! * ST ATI ON* , 2X, 7F10, 4) ) 

750 FORMAT!*! HUB AXIAL FORCE CCEFF. I1000*CFZ) HARMONIC ANALYSIS* 
!/F6.2,*CFZ3*, 10F7. 3/C lOX, 10F7.3I ) 

751 FORMAT!* 0.45R STRESS FOURIER COEFFIC lENTS! 10000*SIGMA/EI • 

1/F6. 2,*SIG • ,10F7.3/!IOX,IOF7.3)) 

752 FORMAT!* 0. 70R STRESS FOURIER C06FFI C I ENTS! 10000*SI GMA/ E) * 
1/F6,Z,*SI37* ,10F7.3/ !10X, 10F7.3) ) 

753 FORMAT! • FORCE COEFF, ! !000*CFZ) VERSUS BLADE AZIMUTHAL ANGLE BY I! 
ICREMENT OF 10 DEG.*/!* CFZ(PSI) *,10F7.3)) 

754 FORMAT!* 0.45R SIGMA VERSUS AZIMUTHAL ANGLE BY INCREMENT OF 10 DEl 
!.•/!• SIG0.45R •♦10F7.3)) 

800 FORMATfSO! ••* )/* RUN NUMBER**, F6.2, 5X, • STRESS PROGRAM OUTPUT CARO 

!*) 

REA0C5.701) ! ( EC! I,UI,J*!.3),I*!,NBE) 
iiRirE(6,701l ( (ECU, J),Js!,3), Is!,N3E } 

READ! 5, 7021 ! f EK f I , J ) , J*1 ,5 ) ,1*1 ,NBE) 

MRITE16,702) I (€<! I,J),Js!.5),I«l,NBEI 
READ! 5, 703 1 IIEM f 1 ,U) , J*!,3 ), I«1,N8E) 

WRITE!6,703) ! ! E Ml I , J) , J*l, 3), I«1 , NBE ) 

R6A0!5,704) ( ( E SI I , U) . J* 1, 31 , 1 «!, nBE ) 

HRITE!6,704I I (ES II , J I , J*1 • 3 1, 1*1, NBE ) 

N04sN8E+! 

N07»NBE*-2 
BRAD * BERCNBE) 

R04»0.45*8ER!NBEI 

R07s0.74<BER!N8E) 

BEL(!l = BeRll )-R0 
00 6 1 = 2, N8E 

6 BEL!I) = BER! I)-8ER! I-l) 

1 WRITE!6,602) 
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REftD«5,705» 0PTl»ITCI}«l»lt9),Q?r2,f T(IS,{slO,l8} 

WRITE<6«&05)OPri* (r«n«I»l«9|«OPT2«(TlI |,I»10* 16) 

1000 REROl 5, 7051 T 
WRITEI6,605> T 
IF( m ) .EO«ENO) RETURN 
IF( T< ll.ME.OATIN) GO TO 1000 

REftDC5»708) RUN, RPM ,NPTS, NPOST,NHARM, FACTOR, 02* I R*I ) =2, NPOST) , 
IRAQI US 

WRITE 1 6, 7061 RUN, RPM .MPTStHPOST, NHARH, FACTOR *02 ,( R<I ) «I «2,NP0ST» 
I, RADIUS 

OMEGA « RPM*3,14159/30, 

O«2s0MEGA**2 


C FACTOR = 0.9625«0MEGA«*2*BER(NBE )*«3 (0,9625 « PI*RH3/4) 

C MULTIPLY FACTOR TO OBTAIN Z/R IN Z 

FACTOR = 100,#FACT0R 
QR » 0.9625*0M24BRAD«'*3^253. 

WRITE I 7, 600) RUN 
NIHARH*NHAR»+1 


R(l )<:R0 


NP0ST1=NP0ST+1 
RINPOSTl ) = RAOIUS 


00 180 I=2,NP0ST1 
180 R( I ) = R( I)*BER(NBE) 


fWARMl»24NHARM+l 
CALL 62C( W7,EK,NBE,5,3) 

CALL 82C(W5, EC,VBE, 3,2) 

CALL LICQM( W6 ,W7 ,1. ,W5,OM2, NBE, NBE I 

CALL B2C(M7,EM,NBE,3,2I 

Mim-0. 


NR=-1 

C 

1002 00 2 NHA»1,NHARM1 

REA0(5,706) lol , ( ril ( I ) , I =2 , NPOSTl ) 

WRITE(6,606) 1 0 1, ( W 1(1 ) , I =2,NP3ST1) 

IF<NHA.NE.l) G3 TO 220 

K1 = I 

CALL MOMEMT(02,l,R,RO,BRAO,0,,Wl,Kl,NPOSTl,KO) 

CT = 0.02*FACTOR*O2(ll/(0R*BRA0) 

CMT =0.01*FACT0R*02I 2) /( QR4BRA0<'»2) 

220 Kl:*l 

DO 3 1*1, NBE 

y2=aER(ii 

Yl*y2-8ELII) 

CALL M0MENTI02,1 ,R,Yl, Y2,Yl,Ml,Kl,NPQSTl ,K0) 

W4(I )*02(2)/BEL( I) ! 

W3m*02ll)*W4<I ) ' I 

IF( UEO.I) GO TO 3 ! 

U2(1-1)*(W3( I)>W4( I-l))«FACrOR 
3 CONTINUE 

W2C NBE )=W4( NBE)*FAC TOR 
N3M*NHA/2 

IF(NOH.EO»NR) GO TO 10 j 

0M2N2*-N0H*0M2 i 

CALL LICOM(H5,W6,l.,W?,OM2N2*NBE,NBE ) ! 

C ELIMINATE FLAPPING ANGLE UNOETERM INAT 13 N AT ROTOR FRE3JENCY 


IF(NOH,NE.ll GO TO 210 

. NBE2S * NBE2-NBE 


1 

DO 200 1*1, NBE 


( 

i 

W5(I*NBE) * 0« 


\ 

200 W5(I^NBE2S) * 0. 


( 

W5INBE2) = 1, 


(f i 
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CALL HINV(»5»N8E«OEr,lWl*!M2» 
WS(MBE2I = 0. 

GO TO 10 

210 CALL NINV<U5*NBE,0ET»IW1»IW2I 

10 CONTINUE 

CALL PROOIW5»NBE«N9E«H2«1*M4) 
DO 11 I = i,NBE 

11 SIGHAf I«NHA)sW4(n/BERCN8E} 

2 NR»NOH 
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WRI TEI&f603r CTtCNT 
00 77 I-l.NBE 

77 HI fn«BERf ll/BERfNBE ) 

N=1 

NHAaO 

WRITEI6.B07I 

WRITE(6,599l * WU I> •! =1 . NBE » 

IFIOPTl.ME* YES) GO TO 78 
WRITE(7,607I 

WRITE! 7,5991 ( W1 ( I) , I >1 , NBE ) 

78 CONTINUE 

DO 50 K=l, NIHARM 
Kl=<-1 
00 51 J»1,N 
fWA*NHA+l 

IFIOPTl.NE, YES) GO TO 51 

WRITE (7, 600) 102 (J) , K1 , < SIG WA( I, NHA) , I=1,N BE) 

51 WRITE (6,600) I 02( J) , K 1, ( SIGMA! I ,NHA) ,I=1,NBE) 

50 N=2 

CALL FSUM3I8,NHARMI ,F2R,SIGMA,NBE) 

WRITE(6,604) 

WRITE(6,599) ( Wl ( II , 1 = 1 , N BE I 

IFOPT2.NE. YES) GO TO 49 
WRITE(7,604I 

WRITE(7,599) ( Wl ( I ) , I =l , NBE I 
49 CONTINUE 
PSI»0. 

DO 52 1=1,8 

WRITE (6, 601) PSI ,IF2R( I ,JI , J=l ,NBE) 

IFnPT2.NE.YES) GO TO 52 
WRITE(7,601) PSI ,(F2Rn »J), J=l,NBE) 

52 PSI=PSI*45. 

CZl = I BB+CF40M2 > / BEL ( I I + 88/ B EL( 2 ) 

CZ2=-8B/8EL(2) 

OR * 100, *0R 

DO 60 NHA=l,NHARNIl 

Nl=NHA/2 

60 FCIl, NHA) = ( CZl*S IGNA( l,NiA)4-CZ24$IGHA( 2*NHA| ) A ( li- (.l ) A*N1 ) /OR 
CALL 62C(W5,ES,NBE,3,2I 
104=1 
107=1 

00 15 NHA»l ,NHAR«l 
C STORE Z/R IN «4 

DO 16 1=1, NBE 
16 H4(I)=SIGMA(I,NHAI/100. 

COMPUTE THE STRESSES IN W2 
CALL PR00(W5,N8E,NBE,W4,1 ,W2) 

11=104 

CALL IPOLl(W2,BER,VAL,R04,NBE,I I, I 04) 

FCC2, NHA)=VAL 
11=107 


C 
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CALL IPOLUt^2 «BES«VAL«R07»N6E* !1» 107} 

FCi3«NHAI=VAL 
0017 1=1, N8E 
17 SI&NAn,NKAI=W2m 
15 COMTINUE 
N=1 
NH4*0 
X0A*0,A5 
X07»0*7 
WRITE16,6101 

WR!TE{6,599) (W1 (II ,Isl,Me£},X04, X07 

00 70K»1,^IHARN 

<1=X-1 

00 71 J=i,N 

NHA=NHA+1 

71 WRITE(6,620I 102 ( J 1 , Kl, (S IGHA( I, NHA I , 1= 1, NBEI , FC( 2,MHA > ,FCI 3,NXAI 
70 N=2 

CALL FSUM3i 8, NUARHl, F2R, SIGMA ,NBE» 

WRirE(6,611 I 
WRITE (7 *61 1 1 

WRITEC6,599J ( Wl( II , I=1,NBE I 
WRITE(7,599J IWI (I | , I = l , M BE I 
PSI=0, 

DO 73 1=1,8 
00 74 I I«6, 7 

74 WR!TE(Il,60U PSI , ( F2RU , J) , J=1 ,MB£) 

73 PSI^PSn-AS. 

00 100 11=6,7 

WRITEUI , 75 01 R0N,( FC( l,NHA» ,NHA*1 ,NHARM1). 

WRITE* 11,751 I RU^^, (F C( 2, MHA ) ,MHA= 1, NHARM 1 I 

100 WRITEUI, 7521 RUN , I FC f 3 , NHA 1 , NHA*1 , '<H ARMl I 
call FSUH3( 36,NHARM1,F2R1,FC,3) 

00 lOl 11=6,7 

WRITEI II ,7531 (F2R1 (1,11,1=1,36) 

101 WRITEUI, 7541 ( F 2R 1C 1 , 2 1 , I =1,361 
DO 102 1=1,2 

REA0I5, 705) T 
WRITE(6,705) T 

102 HRITEI7,705) T 
GO TO 1 

END 



onr> ry nooononoooooooooonooonooooononnoononooooo 


GIRAYIONS DORAKD 



DH-2011 C 
4o 


E 1 


SUBROUTINE ftTElS 
PURPOSE 

COMPUTE THE EIGENVALUES 3F A REAL ALMOST TRIANGULAR MATRIX 
USAGE 

CALL ATEIGIM. Af RR,RItIANA*IAJ 

OESCRIPTION OF THE PARAMEtERS 
M ORDER OF THE MATRIX 

A THE INPUT MATRIX* M 8V M 

RR VECTOR CONTAINING THE REAL PARTS OF THE EIGENVALUES 

ON RETURN 

RI VECTOR CONTAINING THE IMAGINARY PARTS OF THE EIGEN- 

VALUES ON RETURN 

I ANA VECTOR WHOSE Dlf€NS10N MUST BE GREATER THAN OR EQUAL 
TO M* CONTAINING ON RETURN INDICATIONS ABOUT THE WAV 
THE EIGENVALUES APPEARED ISEE MATH, OESCRIPTIONJ 
lA SUE OF THE FIRST DIMENSION ASSIGNED TO THE ARRAY A 

IN THE CALLING PROGRAM WHEN THE MATRIX IS IN DOUBLE 
SUBSCRIPTED DATA STORAGE MODE, 

IA=M WHEN THE MATRIX IS IN SSP VECTOR STORAGE MODE, 

REMARKS 

THE ORIGINAL MATRIX IS DESTROYED 

THE DIMENSION OF RR AND RI MUST BE GREATER OR EQUAL TO M 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

OR double' ITERATION 
REFERENCES 

J.G.F. FRANCIS - THE OR TRANSFORM ATI ON---THE COMPUTER 
JOURNAL* VOL. A, NO. 3. OCTOBER 1961* VOL. A* NO. A* JANJAR 
1962. J. H. WILKINSON - THE ALGEBRAIC EIGENVALUE PROBLEM - 
CLARENDON PRESS* OXFORD* 1965. 


SUBROUTINE ATE I GI H, A .RR.RI * I ANA *I A I 

DIMENSION Am* RRUlf RI (t I.PRRt 2)*PRI( 2)*I ANAI II 

INTEGER P*P1.0 

67=l.OE-fl 

E6=l.0E-6 

ElO«l.OF-lO 

OELTA»0.5 

MAXIT*30 

INITIALIZATION 

N=M 

20 Nl*N-l 
IN=N1*IA 
NN=IN+N 

IFINl) 30*1300,30 
30 NP*N+l 


C 


/ 
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C ITERATION COUNTER 

C 


IT*0 
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ROOTS OF THE 2N0 ORDER MAIN SUBMATRIX AT THE PREVIOUS 
ITERATION 

00 40 1=1,2 
PRRII 1=0.0 
40 PRIIU=0.0 

LAST TWO SUBOIAGQNAL ELEMENTS AT THE PREVIOUS ITERATION 


PAN=0.0 
PAN1=0.0 

ORIGIN SHIFT 

R»0.0 
S=0.0 

ROOTS OF THE LOWER MAIN 2 BT 2 SJBMATRIX 

N2=Nl-l 
INl=IN-I A 
NNl = INl+N 
NIN=IN+N1 
N1N1=IN1+Nl 
60 T=A< NlNl |-A<NNI 
U=T*T 

V»4.0«AININ »»A(NN1 » 

IFIABSi V)-U*E7I .100,100,65 
65 T=U+V 

IFIABSITI-AMAXK U,A8S( VH4E6I 67,67,68 

67 T=0. 0 

68 U=(A< NlNllfAINN)) /2.0 
V=SORTIABSIT ) 1/2.0 
IF* Tl 140,70,70 

70 IFIU) 80,75,75 
75 RR(Nll=U+V 
RR<NI=U-V 
GO TO 130 
80 RR(NU=U-V 
RRINJ=U<-V 
GO TO 130 

100 IFIT)120,110,110 
110 RRfNll=A(NlNl) 

RRINI«A(NNI 
GO TO 130 
120 RRIN1I = A<NNI 
RRINt=A(NlNl I 
130 RHN»=0.0 
RKNl 1=0.0 
GO TO 160 
140 RRCNll»U 
RRIN|=U 
RMNl l = V 
RI<NI=-V 

160 IFIN21l2aO, 1280,180 
C TESTS OF CONVERGENCE 


1 
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180 

RM 0 D=RR«NU*RR(Ml >♦R^Nll♦R HNl ) 

EPS=E 10 »S 0 RTIRM 30 ) 

IFf A 8 S(A(MtN 2 I l-EPS 11280 , 1280.240 
240 1 FCABS(A(MNI n-El 04 ‘ABS(A<NNII) 1300 , 1300 , 250 
250 IF<ABSiPANl-A«NlN 2 l l-ABSJAINlN 2 n 4 E 6 ) 1240 , 1240,260 
2 60 IF(ABS( PAN-A<NNl I 8 -ABS t Al NMl I »*E 6 1 1240 , 1240,300 
300 IF< IT-MAXm 320 , 1240,1240 


E 1 


COMPUTE THE SHIFT 


320 J =1 

00 360 1 * 1,2 
K=NP-I 

IF<ABS(RRIXJ-*PR'U U » ♦ASSI RI f Kl -PRII HI -DELTA*! ABSI RR I Kl I 
I ♦ABSCRKKim 340 , 360,360 
340 J=J+I 
360 CONTINUE 

GO TO ( 440 , 460 , 460 , 480 ) ,J 
440 R= 0.0 
S= 0,0 
GO TO 500 
460 J=N+ 2 -J 

R*RR(J)*RRU 1 
S*RR( J)+RR< J» 

GO TO 500 

480 R=RR(N)*RR(Nl )-RI(N)*RnNl» 

S=RR 4 N)*-RR(N 1 ) 


SAVE THE LAST TWO SU 80 I AGONAL TERMS AND THE ROOTS OF THE 
SUBMATRIX BEFORE ITERATION 


500 PAN=A(NN1) 

PAN1=A(N1N2I 
DO 520 1*1,2 
K*NP-I 

PRR(II*RR(XI 
520 PRI(n*RHK» 

SEARCH FOR A PARTITION OF THE MATRIX, DEFINED BV P AND Q 
P*N2 

IPI = N1N2 
00 580 J=2,N2 
IPI*IPI-IA-1 

IFIA8S(A(IPn)-ePS» 600,600,530 
530 IPIP*IPI+I4 

IPIP2aIPIP + IA 

0=A|IPIPJ*( A« IPI PI-SI *4(1 PI P2)*A IIP IP*U*R 
IFI 0)540, 560, 540 

540 IFIABSIAC I PI ) *Af IPIP+l ) ) *I ABSI 41 IP IP I+AUPI P2*l )-SI ^ABSI All PI P2*2 
1 II -ABSID)*EPSI 620,620,560 
560 P*N1-J 
580 CONTI MJE 
600 0=P 

GO TD 680 
62 0 Pl»P-l 
0=P 

DO 660 1*1, PI - 

IPI=IPI-IA-t 
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IPIABSI A( I 

Pll l•EPS)680,680,660 


66 0 Q« 0-1 

OR OTUBLE ITERATIQM 

680 II*«P“ 1 I*I 4 +P 
00 l?a 0 I*PtNl 
IIlsII'IA 

IFtl-PITZOf 700»720 

700 ipr»m-i 

IPIP»IIP+l 

INITIALI 74 TI 0 N OP THE TRANSFORMATION 

Gi=A(in*i4n n-si*AniP)*AUpn*R 
G 2 « 4 ( IPH*( 4 UPIP) + 4 n Il-SI ’ 
G3»Anpii*AnpiP^i ) 

4 UPI«-ll« 0 .Q 
GO TO 780 
720 Gl=AniU 
G 2 = 4 ( IIl + ll 
IFn-N 2 l 760 f 740 . 760 
740 G 3 *AI III +2 I 
GO TO 780 
760 G 3 = 0,0 

780 C4P*S0RTIGl*GUG2*G2 + G3*G3f 
IF! GAP) 800, 860,800 
600 IF<GI 1820,840.840 
820 CAP-CAP 
840 TsGl + CAP 
PSIl»G2/T 
PS12=G3/T 

ALPHAS2.0/I I .0>PS Il*PSl l*PSI2*PSI 2) 

GO TO 880 
860 ALPHA*2.0 
PSII=0.0 
PSI2=0.0 

880 IF( 1 - 0 ) 900 . 960,900 
900 IFU-P) 920 , 940.920 
920 AUI l)»-CAP 
GO TO 960 

940 Anin=-Anii) 

ROW OPERATION 

960 IJ*II 

00 1040 JsI.N 
T*PSIl*An J+ll 
IFI l-Nl > 980 , 1000. 1000 
980 IP 2 J=IJ +2 

T*T+PSI 2 *A 1 IP 2 JI 
lOOO £T 4 * 4 LPHA 4 <T*AIIJJ) 

All J) *AI UI-ETA 
AIIJ* 1 I*AIIJ*U-PS 114 ETA 
IFII-Nl ) 1020 . 1040, 1040 
102 0 A( IP 2 J)=AUP 2 J)-PSI 2 *ETA 
1040 Ii»I J+I A ' 


COLUMN OPERATION 
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!F( I‘>NU1080, 1060*1060 

1060 

GO TO 1100 
1080 K*H-2 
IlOO IP=IIP-I 

DO 1180 J»0»K 
JlPalP+J 
JIaJIP-I4 
T«PSU*AUIPI 
IFI I-Niui20tll%0*ll40 
1120 JIPZxJlP^lA 

T=T+PSI2*AUIP2I 
llAO ETAaALPHA*IT«'A«Jin 
A< JI MA( JD-ETA 
A( JIP>*AU1 P*“ETA#PSIl 
IFI I«N1 1 1160* 1180* 1180 
1160 Ai J!P2»*AUIP2I-ETA*PSI2 
1180 CONTINUE 

IFU>N211200« 1220. 1220 
12 00 JI=II+3 

JIP=JI4^IA 

JIP2*JIP^I4 

ETA=ALPHA*PSI2*A( JIP2) 

AUH=-€TA 
AC JIPI=-ETA*PSI1 
AC JIP2)=Af JIP2»-ETA*PSI2 
1220 II=IIP+1 
ITalT+l 
GO TO 60 

END OF ITERATION 

12A0 IFC ABSI A! NNl 1 l-ABS ( A(N1N2J> » 1300* 1280* 1 280 

TWO EIGENVALUES HAVE BEEN FOUND 

1280 lANACNiaO 
IANACNU=2 
N*N2 

IF(N2I1A00.1400,20 

ONE EIGENVALUE HAS BEEN FOINO 

1800 RRINI*A(NNI 
RI INI *0.0 
lANACN)*! 

IFfNU14OO*lAOO*l320 
1320 N*N1 

GO TO 20 
1600 RETURN 
END 
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SUBROUTINE 62C(R«B«N«N«H1I 

CE S«P. PERNET OE PASSER D»JNE MATRICE BANOE A UN£ 
MATRICE COLONNE. 

N*NOMBRE OE L IGNES OE LA MATRICE BAND£« : 

M=NOMBRE DE OIAGONALES DANS LA BAMDE. 

M1»RANG OE LA OIAGONALE PRINCIPALE, 

A-MATRtCE COLONNE OE DIM.«N««f2 (SORTIE) 

8 = MATRICE 8AN0E DE OIH-sN^M (ENTREE) 

DIMENSION A( U«B(N.M) 

N2=N**2 
00 I 1=1. N2 

1 A(1I = 0, 

Kl=l 

00 2 [=1»N 

Jl = HAXO<Nl-m«l ) 

IFd.GT.HU Kl=Kl+l 
K=<1 

DO 2 J = J1.M 
IK=(K-1 l#N+I 
1F( IK.GT.N2) GO TO 2 
A( IK)=B( I. J I 

2 K=K+1 
RETURN 
END 




9 ( JUNE 70 ) • OS/ 3 60 FORTRAN H' 


CC'NPILER OPTIONS 


1001 

1 

1002 

2 


TIONS *- MAM£= maINiOPT = 00,LIM£CNT=5 
SOURCE, E BCD I C,NOL 1ST tOcCKtLO 

SUBROUTINE CLEAR(N) • , : 

INITIALISATION DBS CB A 0 ST N=1 
COMMON /BLADE/Ce ( 1 AO) ^ - - :.:,v 

CCMMC:N/TRANS/CT{281 ) ' — ■ 

GO TO U001,1002),N ^ : v.- ^ 

00 1 1=1,140 -- --- - 

CBm=0. . .. ^ ---T;" 

RETURN ■ ' ' ‘ " ” 

00 .2 J = 1,2 81 > 

CT( I )=0. • ■ • ' 

return'. - V -V- 

END 


6, SIZE= 
AD, MAP, 


OOOOK, 

NOEOIT. 


OU DES CT A 0 


!ID,N0XR£F 


N=2; = 





SUBROUTINE EIGMOOIE)^* EC* EKf 0ME6&»0MU8, ')8E » 

COHMON^WORK/WH 2 0> * W2( 201 * W3<20l tW^IVOOi *W5<V00I ! 

0II4ENSIQN EH( I S*ECC U«EK{ n«IWl(20}*IW2(20» ! 

DOUBLE PRECISION OEF«ODC2»,OM i 

OATR OEFNOD/'FLRPPING* »*BENOINS • / 

600 FORN ATI *l BLADE E IGENMOOES AT ONEGA **#F7«2** FOR'tHt* ELEMENTS*//' 

I* ROTOR FREOUENCT* *F7.2»* CPS»» SX**R0T3R RPM*tF7#l//" ♦♦RESONANCE 
2 FREOUENCT OF UNIFORM BLADE IFUBI • ,F7*2** CPS*/// 1 

33X*» MODE* *9X,*E16ENV ALOES*. l2X**EIGENFReOUENCIES*,3X. *F/FUB*/ ' 

4 17X. •(RA0/S)*«2« /14X.*REAL**8X.* IMAG*/» 

601 FORMATIIX.AS* 12* IP2E12.3* 8X.0PF7, 2*8X»F6* 2» I 

FUB = OMUB/6.2832 

FROT = 0M6GA/6.2 832 
RPM = 60*4FR0T 
CALL B2C(H4«EK*NBE*5*3I 
CALL 82C< W5*EC.N8E* 3, 2J 
X2=+OM6GA*^2 

CALL LICOHf U4*H4 «N5*X2, NBE«N8E) 

CALL 82C(K5.EM*N8E* 3*21 

CALL MINVrW5.NBE.DELT* IWl. IW21 

CALL PROOUW5*Wi^*Wl*N6E*NBE) 

WRIT EC 6* 6001 OME 6A , NBE , FROT ,RPN *FU8 
CALL HSBGINBE*W4*NBE» 

CALL ATEIG(NBE*W4»WI*W2»I Wl.NBEt 
DM = DEFHODI I i 
1 001 DO 1 I*1*NBE 
MODE = I-l 
K = NBE - MODE 

FREO « S 0 RTrA 8 S(WllK)l i/6.2832 
FREOR = FRE0/FU3 

WRITE(6,601 I DN.MO3E*Ml(<I.W2l<J,FRE3.FRE0R 
1 DM = DEFM0DI2I 
RETURN 
END 
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SUBROUT INE FHARM3(NPTS.NHARf^tNPOST tWTYPE) 

COMMON FOi 20 J ♦FN41 2.20,21 * F« 24, 20, 2) 

COMPLEX EJO.EJND.CN 

N1 = NPTS + I 

PI=3. 141592 

PIDX=PI4*2/NPTS 

EJO = CMPLX<0.,-2,*PI/FL0ATINPTSn 
EJD = CEXPIEJDI 

1010 DOlO I=1,NP0ST 

Fom=o, 

1001 001 K*1,NPTS 

1 FO(n*FO<n+F«K,I,NTYPEI 

10 F0CII = F0<n/NPTS 
EJN0=U.,0.1 

1002 002 N=1,NHAP.M 
EJND=EJNO*EJD 

COR = U, - REAL(EJN0n/4 PIDXX«FLGAT(N#<‘2J » 

1011 DOll I=1,NP0ST 
CN=(0.,0.l 

1003 003 K=1 ,NPTS 

3 CN=CN*EJNO+F(Nl-K, ItNTVPE) 

CN = CQR*CN 
FNIN.I,1)=REAUCN1 

11 FN(N, I, 2> = -AIMAG(CN ) 

2 CONTINUE 
RETURN 
END 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


SUBROUTIN6 HSBG 
PURPOSE 

TO REDUCE A REAL MATRIX INTO JPPER ALMOST TRIANGULAR FORM 
USAGE 

CALL HSBG{N«A«IA) 

DESCRIPTION OF THE PARAMETERS 
N ORDER OF THE MATRI X 

A THE INPUT MATRIX, N BY N 

lA SIZE OF THE FIRST DIMENSION ASSIGNED TO THE ARRAY 

A IN THE CALLING PROGRAM WHEN THE MATRIX IS IN 
DOUBLE SUBSCRIPTED DATA STORAGE MODE, IA=N WHEN 
THE MATRIX IS IN SSP VECTOR STORAGE MODE. 

REMARKS 

THE HESSENBERG FORM REPLACES THE ORIGINAL MATRIX IN THE 
ARRAY A. 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

SIMILARITY TRANSFORMATIONS USING ELEMENTARY ELIMINATION 
MATRICES, WITH PARTIAL PIVOTING. 

REFERENCES 

J.H. WILKINSON - THE ALGEBRAIC EIGENVALUE PROBLEM - 
CLARENDON PRESS, OXFORD, 1965, , 


SUBROUTINE HSBG(N,A,IA» 

DIMENSION AllJ 
DOUBLE PRECISION S 
L'N 

NIA=L*IA 
LI AaNlA-IA 

L IS THE ROW INOEX OF THE ELIMINATION 

20 IF(L-3I 360,40,40 
40 LIA=LIA-IA 
L1=L-1 
L2»L1-1 

SEARCH FOR THE PIVOTAL ELEMENT IN THE LTH ROW 

ISUB=LIA+L 
IP I V= I SUB- 1 A 
PIV=ABS(A( IPIV II 
IFIL-31 90,90,50 
50 MsIPIV-IA 

DO 80 I»L,M,IA 
T-ABS(Af II I 
IF<T-PIVI 80,80,60 
60 IPIV*I 



ono nno ono r> o o o n o non 


PiVsT 

80 CONTINUe 

90 IF«PIV» 100,320*100 
100 IF(PIV<>ABS(AI ISJB m 1B0*180»120 

INTERCHANGE THE COLUMNS 

12 0 M=IPIV-L 

03 140 l*l*L 
J=M+I 
T=A|J» 
tCsLIA+l 
A<J)sA(Kl 
140 A<K| = T 

INTERCHANGE THE ROWS 
M«L2-M/IA 

00 160 I-L1«NIA,IA 
T=AI 1 ) 

J=I-M 
A(TMA( J} 

160 AIJ»=T 

TERMS OF THE ELEMENTARY TRANSFORMATION 

180 00 200 I=L,LIA»IA 
200 A( n-A( n/AI I SU3I 

1 

RIGHT TRANSFORMATION 
J=-I A 

00 240 I»ltL2 

J=J+IA 

LJ=L+J 

00 220 K=l*Ll 
<J=<+J 
Kt=K+LI A 

220 AIICJI=AIKJI-A<LJI*A(KLI 
240 CONTINUE 

LEFT TRANSFORMATION 

K=-IA 

00 300 1*1, N 

K*K+IA 

LK*K+L1 

S*A1LKI 

LJ=L-IA 

00 280 J*l,L2 

JK=iC+J 

LJ*LJ+IA 

280 S*S+AILJI*AI JKI* l.OOO 
300 A(LKI=S 

SET THE LOWER PART OF THE MATRIX TO ZERO 

DO 310 I*L.LIA,tA 
310 A(I}*0.0 
320 L*Ll 

GO TO 20 
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SUBROUTIME IPOLl < Y« X« VAl ,ARG«N« 1 1,1 2) 

SOUS PROG. D» INTERPOLATION L INEAIRE 0»0NE FONCTIDNY OONYEt 
PAR Yin ET X(N .ARGsABSCISSE 0» INTERPOL. ET VAL LA YALEOR 
CORRES PQNOANTE. 

N=NOMBRE OE POINTS DEFINISSANT LA FONCTION. 

Il»VALEJR 06 DEPART OE LA VARIATION OMNOICE. 

I2»IND!CE OE LA PREMIERE VALEllR SJP. 00 EGALE A ARG. 
oiHENsiQN xin ,vm 

00 I 1=11, N 
12 = 1 

IFiXi n.GE.ARGl GO TO 2 

1 CONTINUE 
GO TO 3 

2 IFi I2.E0.n GO TO 3 

VAL=IYU2»»< ARG-XI I2-ll»«-YC 12-11*1 XU 21-ARGll /I XlI2l-xn2-U I 
RETURN 

3 VAL = Yn2J 
RETURN 
END 



non 



SUBROUTINE LICQH IX, XI . Cl » X2 9 C 2 * N* N 8 

CE S*P. EFFECTUE LA COMBINAISON LINEAIRE DE DEUX 'lATRICES 
C3LONNES XI ET X2 OE 31 M. *N«M» AFFEC TEE S RE SPECTI VEMENT 
DES COEFF* Cl ET C2.LE RESULTAT EST OQNNE EN X* 

DIMENSION XHU,X2IU,XCU 

«N*N*M 

DO 1 I si, MN 

1 XI I 8»C1*XIII l+C2«X2C 18 
RETURN 
END 


i 

I 
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SUBROUTINE MI NV 


PURPOSE 

INVEilT A MATRIX 


USAGE 

CALL MINVI A,N,0,LtM> 

OESCRIPTION OF PARAMETERS 

A -"INPUT MATRIX, OESTROTEO IN COMPUTATION AND REPLACED BT 
RESULTANT INVERSE. 

N - ORDER OF MATRIX A 
0 - RESULTANT DETERMINANT 
t - WORK VECTOR OF LENGTH N 
M - WORK VECTOR OF LENGTH N 

REMARKS 

MATRIX A MUST BE A GENERAL MATRIX 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

THE STANDARD GAUSS- JORDAN METHOD IS USED. THE DETERMINANT 
IS ALSO CALCULATED. A DETERMINANT OF; ZERO INDICATES THAT 
THE MATRIX IS SINGULAR. 


SUBROUTINE Ml NVI A ,N »0 , L, Ml 
DIMENSION Al II ,U II .Ml II 


IF A DOUBLE PRECISION VERSION OF THIS ROUTINE IS DESIRED. THE 
C IN COLUMN I SHOULD BE REMOVED FROM THE DOUBLE PRECISION 
STATEMENT WHICH FOLLOWS. 

DOUBLE PRECISION A. 0, BIGA, HOLD 

THE C MUST ALSO BE REMOVED FROM DOUBLE PRECISION STATEMENTS 
APPEARING IN Ot-iEft ROUTINES USED IN CONJUNCTION WITH THIS 
ROUTINE. 

THE DOUBLE PRECISION VERSION OF THIS SUBROUTINE MUST ALSO 
CONTAIN DOUBLE PRECISION FORTRAN FUNCTIONS. A8S IN STATEMENT 
10 MUST BE CHANGED TO DABS. ! 


SEARCH FOR LARGEST ELEMENT 


0 = 1 . 0 

NK»-N 

DO 80 K=l.N 

MK=N<+N 

LIKI»K 




non oooo ooo oor> 


GIBA7I0NS DOHASB 

HIKI=K 
KH=N<+K 
BIGAsA(KKI 
00 20 J=K,M 
U=M»{ J-1 i 
00 20 I»K,N 

10 IFI ABSfBlGA}- IVBSMUJH) 15t20»20 
15 BlGAsAdJl 
U K |« I 

2 0 CONTINUE 

INTERCHANGE ROWS 

J*L(K| 

IFIJ-KI 55*35*25 
25 KIa<-N 

□0 30 1 = 1. N 
KI»KI+N 
H0L0=-A(KII 
JI*<1-K+J 
A( KI)=A( jn 
30 A( JI I =H0LD 

INTERCHANGE COLJMNS 

35 1=M(KI 

IF(I~K> <>5.A5*38 
38 JP*N*( I-U 
00 40 J*l*N 
JK=N<+J 
JI*JP+J 
HnLD=-A(JKI 
AU<} = A|JII 
40 AUll “HOLD 

DIVIDE COLUMN BY MINUS PIVOT (VALUE OF PIVOT ELEMENT IS 
CONTAINED IN BIGA) 

45 IF(BIGA) 48,46,48 

46 0=0,0 
RETURN 

48 03 55 1=1, N 

IFII-KI 50,55,50 
50 IK = NK-»-I 

A(IK}=A(IK|/<-BlGAI 
55 CONTINUE 

REDUCE MATRIX 

DO 65 1=1, N 
IK=NK+I 
IJ»I-N 
00 65 J=l,N 
IJ=IJ+N 

IF(I-K» 60,65,60 
60 IFIJ-KI 62,65,62 
62 XJ=IJ-I+K 

AdJ l = A( IK)«A(KJ IJ t 
65 CONTINUE 
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DIVIDE ROW BV PIVOT 


KJs<-N 
DO 75 Jsl,N 

IFIJ-KJ 70t75.70 
70 A( KJ) «A(KJI /BIG4 * 

75 C 0 » 4 TINUF 

PRODUCT OF PIVOTS 

0*0*BIGA 

REPLACE PIVOT BY RECIPROCAL 

AIKICI® 1 , 0 /BIGA 
80 CONTI NUE 

FINAL R 3 M AND C 3 LU'^N INTERC-IAN 56 

K*N 

100 K=IK-U 

IFIKl 150 ,l 50 fl 05 
105 I»LKJ 

IFU«KI 120 , 120,108 
108 J 0 *N*IK- 1 J 
JR*N*( I-l ) 

00 110 J*1 ,N 
JK = J04-J 
H 0 L 0 *A(JK) 

A<JO=-A(jn 
110 AUn =HC «.0 
120 J=M<K) 

IF(J-K> 100 , 100,125 
125 KI*K-N 

00 130 1 = 1 , N 
KI=KI + N 
H 0 LD=A(KI) 

JI = <I-K+J 
A(Kll=~AUn 
130 Ai JIl =H 0 L 3 
GO TO 100 
150 RETURN 
ENO 
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SUBROUTINE MQNENTt F XN« Nf X »)( 1«X2 «XOR«F «K1 , K2t KO ) 

C TO COMPUTE THE MOMENT UP TO THE NT^ ORDER OF A MONOT3NICALL Y TABULATED 
C FUNCTION F.IN THE INTERVAL f Xl« X2) «LE« <X (K1 ) » X(K2 )1 WITH RESPECT TO PO 
01 NENSION FCUtFXNd i«X(ll I >SoNT XOR 

lOOl 00 I K*Kl,K2 

C SEARCH FOR FIRST VALUE OF X LESS THAN XI 
IFIXIKI.GT.Xll GO TO 2 
n»K 
12*11 

C SEARCH FOR FIRST VALUE OF X LESS THAN X2 

2 1F<XIK).LT.X2I GO TO I 
I2=K-l 

GO TO lOOT 
; I CONTINUE 

1003 00 3 M*l,NM 

3 FXNiMlaO, 

1004 00 4 I *11.12 
J=I + l 
X1*X(I1 
XJ*X<J» 

IFIUEO.m XI=X1 
IFIl,E0,I2l XJ*X2 
A=Fin 

C*|F<JI-AI/<X(JI-X( m 
Zl«XI~XOR 
Z2=XJ-X0R 
A=A+C*(XQR-X( in 
2004 00 4 M = l.NM 

Z12=Z14I XI-XORI 
Z22*Z2*|XJ-XORI - 

FXN( MI=FXN( M»+A*( Z2-ZU /M+C4( Z2 2-Z12 I / IM*l * 

Z1*Z12 

4 Z2*Z22 
KO*I 1 
M = I2 
RETURN 
ENO: 
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SUBROUTINE PROO( A«L«N« A6I 

aiNENSION AUl «BC1I «ABm 

00 I I»l,M 

iM » n>i)*N 

00 2 J-1«L 

j I = j ♦ n - u *L 

ABUn - 0. 

00 2 K«1,N 
IU = <*!M 
JKaj4-(K-Un 

2 AB(JU » AB4JU 4- A(JK)«3(Kn 

I CONTINUE 
RETURN 
END 
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SUBROUTINE PR001«A.B,X*N,M) 

COMPUTE THE PRODUCT A* B OF Tk«3 N»N AND N*M MATRICES 

STORED COLUMNWISE AND STORE IT IN B 

DIMENSION AUI*Bfl)tXm 

DO I 1 = 1, M 

IM»U-n*M 

00 2 J=l,N 

XU 1*0, 

00 2 K=l,N 
KI=K+IM 

2 XUI=XU)4'AUfa«Bf XU 
00 1 K=I,N 
KI=K+IM 

1 B{KU = X|XI 
RETURN 
END 
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SUBROUTINE STRftHI N96*NR* I S*ROt StlO #HE ,R S*EMt ES* EC tEiU 
COMMON/WGRK/Kl (23) *W2<20» *W3I20»t 0(^03)* C(400) 

CONN0N/BtADE/BEL(2O) *BER(20)*BECN(20) tBETM(20) *BEK(20I# ; 

I A(?0lf 8(201 

OIMENSION ESCNBE *31 •EM(NBE»3)*€CINB6» 3 )» EK(NBE»5 )t HE< 1 »*RS( Ilf 
1 8110(4} ' 

600 F0RMAT(/IX*4A4* 13* • EL EMENTSS 3X* •CF=« * IPE 10* 3 * 3X* *BB-* *E 10* 31 | 

700 FORMAK ♦ IBLAOE IDENTIFICATION • , VA4*2X**NBE*** 13* 3X* «CFa« * IPEIO* 3*| 

!• 8B«**E10*3/( • BER**6X*7E10*3} I | 

701 FORMAT (• EC* * 7X# IP3 EIO.3) ! 

702 FORMAT!* EK* * 7X* 1P5E10.3) 

703 FORMAT! • EM* . 7X* IP3E 10*3) 

704 FORMAT!* ES* *7X*IP3E10.3 » 

N8El» N8E-1 

f NBE2 a NBE-2 

EMd.lt = BECM(l)/4. 

BEUNBE+1) a 1. 

U!l ) a 1./9EU 1 I 

RADIAL DISTANCES FROM THE BLADE ELEMENT CENTERS ARE STORED IN Ml 

WKII a BELIU/2. ♦ RO 
IF(N8E,E0.l) GO TO 1203 

1202 DO 202 1=2. N8E 

Wim = BER!I-U ♦ BELm/2. 

202 U( n a l*/3EL(I) ♦ l*/BELU-l} 

COMPUTE AfB AND C. (STORE THE HI^RI IN W2' AND W3» 

W3! NBE+1 )=0* 

1203 DO 203 I=l*N8£ . 

H2(U a BETM! n*8ER( I ) 

W3 (U a BECM(I}*Wl (I I 
A( I) a W2m + W3(n /2. 

203 Wl n ) a BEX! I»*U< I > 

Ml CONTAINS NOW THE PRODUCTS BEK^’J 

WUNBE-ni a 0* 

CINBEJ a 0* 

IF!NBE.EO*ll GO TD 1206 

1204 DO 204 1=2. N8E 
X=N8E-I-»-2 

204 C(<-1 ) a C!K» ♦ M2!X} ♦ M3(K} 

COMPUTE EM a MASS MATRIX 

EC a CENTRIFICAL MATRIX 

EK * RIGIDITY MATRIX. 

1205 DO 205 I a 2.N9E 

ES(I-l*ll a l.7(BELin*BELC I-IM 
ES! 1-1.2) a -U!II /BEUn 
ES(I-1.3I a l./BEL(n**2 
EKd.U a BEK(nFES!{~l«U 
EX II. 2) a -(wim+wi( i+in/BELin 
ECn.U a -(AID^CU n/BEL( I» 

205 £M(l. II a BECMd ) A. 

c FREE HINGE 

ES(l.l) a 0. 



uuu uoo ou uo 



00 226 

ECU. 3t EC4 Ui.li 
E^«II.3} = EM< IH.ll 


22 6 EKU ,61 * EK< 1>1 ,2* 

IF* J^BE*LT.3I GO T0 208 

00 207 I»1,NBE2 

207 EKU.5I « EKI U2 ,1» 

208 CONTINUE 

1206 DO 206 !»ltNBE 
Bf n»W3U«'U/2« 

EKII,3i * BEK4 II /BELtl ♦ «U I»U*UI H-ll ♦ BEKIU2>*ES<I,3) 

ECII.2) * C4IKUIUII *■ A( n/6ELl I S “ BU }/BElf 1*11 
206 EMU, 21 =: EMUfli > EMU, 39 ♦ 8ETM<U 
COMPUTE TOTAL CENTRIFJC4L EFFECT 

CF * Clll ♦ 8ECM<l)#CBERU»-!-B£L{U/2.9 ♦ BETMt 1 9 *BER4l ) 

II = I 
ARG * RO 

R106 a 5000*B6RCNBE I 

COMPUTE THE ELASTICITY MATRIX ES, 

10000#( SIGMA /El = ES*IZ/8RA0l 

1222 DO 222 1=1, NBE 

CALL IP0UIHe,RSfVAL,ARG,NR,Il»I2l 

VAL = VAL*R106 

ESII ,U = ESU ,n*^AL 

ESI 1,21 = ESI I,2l*VAL 

ESI I ,31 = ESU,3 ^♦YAL 

ARG = BERIIl 

222 II » 12 

COMPUTE AXIAL FORCE COEFFICIENT DUE TO BLADE BEN01NG=38 
88=BEKI 2»/8ELIli 

IFIIS.EO.O) RETURN 

WRITE<6,600I BLIO.NBE.CF, 88 

PUNCH THE RESULTS 
2223 DO 223 K=6,IS 

IFI<.EQ.6I GO TO 7701 

MRITEIK,700I BL ID, NB6, CF, 68 ,R0, I SERI 19, 1 = 1 ,NBE I 
7701 WRITE<K,701I I ( EC II , J9 , J = 1 , 39 , 1 =1 ,NBE I 
WRITE (K, 7029 I ( EKII , J 9 , J = 1, 59 ,1 =l ,NBE 9 
WRITE(K,703I (IEMII,J9,J = 1,3I, I*1,N8E I 

223 WRITEIK,706| f ( E SII , J9 , J = 1 , 31 ,I =1 ,N8E 9 
RETURN 

END 



O.^SRCOMPUTED ^TR£SS£3 COURIER COERR/C/BHT^ 


6~ in hactobar^ 

fz_£f_ 

® I (D I ® I (i 


cos'f + <T, + (S', 




0.09 


0.25 


-0.04 

0.14 

0 , 2 a 

-0.36 

-O.OS 

-OM 

-0.66 

6.21 

WSSM 


0 09 

wrnm 

^^9 

-0.2 

-0,24 

0.3 

-O.SG 

-o.3e> 

0.30 

-0.31 

O.H>5 

-0.04 

0 

0.12 

-0.09 

Bl^9 

013 

-0.4G 

-0.62 

0.11 

RRSI 

-1.34 \ 

0.6Z 

-0.12 

0.06 

0.76 

6.29 

0.11 

o.a3 

0.72 

0.05 




















































































0.45 R MEASURED STRESSES Toc/R/ER COaFFfCfENTS 


< 3 " in h<Lc.tbh<ir^ 


G"= (3"o 4- C05 4'* 4- 6^, + (T^ <ix>s 


/ 

S’K + G'<, 2 4^ + 



@ 
















9 . 05 


9 . 0 ^ 


f2.ro 


12.11 


12.12 


12.13 


14.10 


14. If 


14 . 12 


14.13 


16 . 08 


16.09 


1646 


8.25 


7.63 


7.63 


C .369 


0336 


. 0.196 


0. 670 


0.213 


0 . 66(3 


.27 


- 0.006 


0.556 


0.911 0.924 


0 652 0.641 


0.224 


- 0 . 547 \ 0.0866 


- O 722 - 0 . a 25 


- 0.0219 


0.396 


- 0.429 


- 0.333 0.0697 


-0.168 -0.269 


- 0.414 I - 0.631 


- 1.01 


- 0 . 7 ff 


0.0515 


0.556 


0.435 


0.260 


0.115 


2.38 


2.40 


2.24 


2.33 


2.68 


2.92 


1.24 


2 . 


2.69 


3 . 


0.77 S 


. 1.06 



0.718 I- 


0.323 


2.62 


3.69 


0 226 


2.22 


.92 


5.91 


0.006 


0.777 


1.79 


1.3 6 


1.06 


7.6 


2.37 


2.17 


. 1.19 


0.856 


. 1.05 


0566 


0197 


0.633 


0.641 


0.365 


0.466 


0.595 


0 255 


- 0.660 


0.926 


0.165 


2.fd 


2.15 


0J6ZS 


GIRAVIONS 
DO RAND 

























































































































GIRAVIONS 

DORAND 


BLodo. sS<LC,C/‘oO 
inarCrct 



miiLwtff.ui'ii 
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^-PROQ.RA«._.ap,II ON# .CG.N^P.UTED values ^Ofr„BLAD£^,£ -- 

-?,pp QgIa^IqFt lO^N 


BLADE STRUCTURAL DATA 

ELE^ENf .. : : . ..^ Ad, H lk.GZ»i.l7:S£CXiXNEfit7i^^ 


2 _ 1.0.COE_QQ._ ,6,500e-02 „_5.«O.QQE... 01 5.000E-06. 

'7^^ sooei^oo .r' ^»50 qe-Q2.c-^£.oooe 

__A 3..0C0E-00 6. 5,0QE-02.,^^,,5, QOOE 01 -5. P00Ei?.06. 

- -'-I. _. 4 . 000 e 00 . ■ 6, 5 OO^-O 2“ -5, bdOE 01- - -:-5 . OOOE- 06 

6 5.G00E- 00 6.500E-02-- -5.000E 01 5. 0006-06 

‘ 6.0C0E 00 6,5006-02 - -5. OOOE 01 5.0006-06 



RADIUS 6. CODE 00 METRES V--.: : - w. ^ _ 

total mass 2.842E 02 KILOGRAMMES _ 

-AVERAGED SECTIONAL INERTIA ... 5.000E-06 METRES##^ 

UNIFORM BLADE FREQUENCY 6.739E 01 RAD/SEC 
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SZi 2..''R(H]QRlEREaullSE^ 


I^feEQM Nt¥:2iiZLJ!aiEjaRMl&LADE; 


7 77- ” E 7l^ E I G E N V A jUJE I G E N FRE QU.EN Cl E S , E/£U G 


- - - - ... - - 



__ . BE.AL 





FLAPPING 

’~0 

’ , 8 . 6.4 5 E 01 

0,0 

■"■••-— 

1.48. . 

0. 14 - - 

IiI'2”B£NQING 

. 1 

4.576E03 

: . 0.0 

,;,Z . 

.L^^.10.77 

-..a. .-00.-- 

BENDING 

z 

4.780E 04. 

0^0-,.^ 


^_-34.80 

. 3 *24 

^ “ ft'p N n T M A 


'■p f fi! E 0«5 

0-0 ■ 

■ r-iT 

.L:.-.V:’72“;.60--'- ■ 

' 7 ' ' • 6 V 7 7"'^ 

BENDING 

4 

.„6.C82£ 05 

_ 0.0 --- 


. . 124.12--, 

U .57 

■.: ... ._. BENDING 

5 

- 1.415E 06 

0.0 ;■ 

. ...; .•:. ■ -. .;■ .; .:. . .• ‘ 

189.32 

■••■■••- 17.6^^:^^--- 

B.ENDING 

6 

. .2. 83BE 06 

0.0 

„ . . 

. 268.11- . 

-25.00— 

^21. -u '.BENDING 

7 

5.1256 06 

0.0 - 


360.29 


BENDING 

8 

8.5 52E 06 

.. 0,0 



.465.42 

. 43.40 

".' "'.F BENDING 

9 

1.34 IE 07 

0.0. .. 


582.78 

54.34 —.^:"^-- 

BENDING 

10 

1.996E 07 

0.0 


.711.09 

66.30 

. BENDING 

11 

2.84CE 07 

0.0 


: 848 .18 

79.03 

BENDING 

12 

3.E74E 07 

0.0 


990.63 . 

92.37 - 

...-■BENDING 

13 

5.C7CE 07 

0.0 •■■ 


-1133,29 

105.67 — -:.7 : 

bending 

14 

6. 357E. 07 

. .0.0 


1268-9.7 

118.32 

.BENDING 

15 

7.613E 07 

0.0 


1388.68 

'■ 129.48 -. 

BENDING 

16 

8.679E 07 

0,0 


1482.74 

138.25 

'■ ■ . BENniNG 

17 

9- 3966 07 

0.0 


1542.73 

143.84 .._ .. 

♦PROGRAM 

OPTION* EIGEN 

MODES OF 

THE BLADE 

» ■ 

l..r 2'.. 1 -. 
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=at,A QE ;E 1 G Eii.MQ d ESl. AT- OM E Gl: l=~:l 3 


!i;I-lrr:__RQT‘D]i[„ £REd:UENtY^r''^ .K3d 


B I G ENV A LUES,..^.^^^ E I.GE.N £R£Q.U ENC I ES _ . F /FUB . 




REAL. _ . 

---: . - 

_ .I.MAG . . 



- _ -™. :z. '. -.^v:: — ‘“. j:; '. 


'1 r"'-: ^ vj^' ' 7 " 'T' 

... ... ' ; "-“H'”-.”!’ 

FLAPPING 0 

1 . 26,0 E. 

03 

.... 0 . 0 . 

5 . 65 _.. 

... 0 .- 53 _ 

i:zr--: ;z 3 ending 

1 

1 . 220 E 

04 .. 

• 0 . 0 ,,. :.......-..'^;±._:..£. 

- _:: 17.58 i.. .£ . 

- • £_„1 ,.6 4 ' - 

BENDING 

2 ,. 

6 . S 34 E. 

.04 

0.0 

-—. 41.91 

-. 3 . 91 . 

■ r::.ia:'ENDlNG. 

3 

2 . 52 IE 

05 

... 0^0 1 .' 

___ 79 . 92 -:— : 


BENDING 

A 

6 . 841 E 

05 . 

. 0 . 0 ., 

131.64 

- 12 . 2.7 

:.L.-.-l. BENDING 

5 

U 533 B 

06 

■ 0.0 . ■■ ■ 

197.05 - 


..._ . BENDING 

6 

3 . CC 9 E 

06 

0.0 

276 . 09 . 

. 25.74 .. . 

r ... ...BENDING 

7 

5. 3626 

06 

0.0 ......... . 

368.54 

34 . 36 ..:....^-— 

_ . BENDING 

8 

8 . 870 E 

06 . 

. , 0 . 0 . 

473.99 

44.20 

• ■^'■;;; bending 

9 

1 . 382 E 

07 

0.0 

591.74 

55.17 

BENDING 

10 

2 . C 49 E 

07 

0-0 

720.44 

67, 17 . 

bending 

11 

2 . 906 E 

07 

0.0 

857.98 

80.00 

BENDING 

12 

3 . S 55 E 

07 

0.0 

1000.89 

93.32 

■bending 

13 

5 . 167 E 

07 

0.0 

1144.00 

106,67 ■: 

BENDING 

lA 

6 . 469 E 

07 

0.0 

1280.08 

119.36 

BENDING 

15 

7 . 739 E 

07 

0.0 

1400.12 

130.55 : 

bending 

16 

8 . 817 E 

07 

0.0 

1494.46 

139.34 

••' BENDING 

17 

9 . 553 E 

07 

0.0 

1555.58 

145.04 


♦PROGRAM OPTIDN« CONTINUE 

♦PROGRAM OPTION* END OF THE COMPUTATIONS 






















Rl/K 9.C3 

CCKFARISCN BETkEEN 
STPESSES At .A5*R 
C.C 

PSI 1 

0 1 

-1 

-I 

30 1 — 

. - 1 
-1 

60 —“I 

-I 

-1 

90 -^-1 — 

-I 

-I 

12 C 1 — 

-I 

-1 

150 ~I — 
"I 
-1 

180 1 — 

-I 

“1 

210 -“-I — 

-I 

-1 

2 A 0 1 — 

-I 

•I 

2?C -^-I — 

-I 

300 1 — 

-I 

-I 

330 1 — 

-I 

-1 

CCKPUTEC STRESSES 
MEASURED STRESSES 


CCMPiTEO STRESSES!*! IFRCR ERC SIMLATICRS) AKO PEASUREC STRESSES C *1 l« ARCH 71 TESTS! 




IE CC 


O.IAOE 01 


I 

CCRPUTEC 

MEASURED 

--“♦I 

C.130E 

01 

C.137E 

01 


G«130E 

01 

0*1AOE 

Cl 

♦ 

0.128E 

Cl 

0*139E 

01 

♦-—I 

C.123E 

01 

0.132E 

01 

1 

O-llTE 

Cl 

0«119E 

01 

I 

0«111E 

Cl 

0.1C2E 

01 

^ ^ 1 

O.IOCE 

Cl 

0,E59E 

CO 

1 

0.103E 

Cl , 

C.766E 

00 

1 

O.IOIE 

01 

0-777E 

CO 

— -1 

C.IOCE 

Cl 

0.E93E 

CC 

1 

C.9S6E 

CO 

O.IOTE 

Cl 

1 

C,95TE 

CO 

0.122E 

01 


C.S13E 

CO 

0*12@E« 

01 

1 

0-857E 

00 

0«121E 

01 

I 

C.BC3E 

CO 

C.997E 

CC 

J 

C.766E 

00 

O.TieE 

CC 

I 

0.756E 

00 

0*A6AE 

00 

I 

C.17TE 

CC 

0.326E 

00 

1 

C«e2AE 

CO 

0-355E 

CC 

I 

C«B81E 

CO 

0-5A3E 

CC 

1 

C.93AE 

00 

0.825E 

CC 


c,s«ie 

CO 

O.lllE 

01 

1 

C.^IAE 

€0 

0«130E 

01 

4 1 

C,958E 

CO 

0.13AE 

01 

1 

0.92IE 

GO 

0.12SE 

01 

1 

C.897E 

CO 

0-IC7E 

01 

1 

c.eicE 

CO 

C.tTSE 

CC 

I 

C.873E 

CO 

C,743£ 

CG 

I 

0.69 IE 

GO 

O.TCOE 

00 

1 

€•9276 

CC 

G.7C7E 

CC 

— - 1 

C-977E 

CO 

0.885E 

CC 

1 

O.IOAE 

Cl 

G.1C2E 

01 

I 

C.llOE 

Cl 

0.113E 

Cl 

— 1 

c.iieE 

01 

0.122E 

01 

1 

C«122E 

€1 

C.127E 

01 

4 1 

G.127E 

Cl 

0.132E 

01 


ARE CIVICEC BY THEIR PEAR YALtEs 
ARE CIVIDEC BY THEIR PEAK VAttE* 


A.790HECT0BARS 

1,93CHECTCBAPS 


OIBATIOI^ » BH^O1 1 • C 



RUK 5 ,04 

CCfFAPlSCN etTWtEN CCMPtTEC 
STRESSES AT .<i5^R 


STRESSES(*IIFRCf' EK SIRtLATICNSI AKC MASUPEE SIRESSESM) (MARCH 71 TESTS# 


C.C 


C-6SIE CC 


0. 


PSI 

0 


3C 


€C 


^0 


12C 


150 


180 


210 


24C 


270 


3C0 


330 


I 


ir ■ f n -t- 

I 

CCRPUTED 
C.118E Cl 
C.lllE Cl 

MEASURED 
C,ia2E Cl 
0«123E 01 




♦ 4 I 

-1 



.44 1 

C.IUE Cl 
C.116E Cl 
0.116E Cl 

C.125E Cl 
Ool2AE Cl 
C.118E 01 




44 I 




♦ « i 

C.llEE Cl 
C«116E Cl 
C.llEE Cl 

C,1C7E Cl 
C«SA4E CO 
O.eSSE CC 



♦ 

4 I 

*1 


4 

4 I 

C.115E ei 
C.113E Cl 
0.109E €I 

0.852E CO 
C*S51E CO 
0«112E 01 

~z\ 



44 I 

- 1 



4 4 t 

C.IOAE Cl 
C.979E CO 
C.?1«E CC 

C.129E 01 
C«138E 01 
0.1338 Cl 



♦ 

4 I 

-I 


t 

4 I 

C.ESAE CC 
C.EOTE CC 
C,TCCE CO 

C.112E Cl 
0.6C7E CO 
0.5C1E CC 


4 

4 

1 

1 

4 

,4 


C.775E CC 

C.3C<^E CO 

1 




C,752E CO 

0.2E8E CC 

^ I 

4 

4 

I 

C.E22E CO 

C.ACCE CC 

^ 1 


4 4 

1 

C.858E CC 

0.7 see CC 

* ^ ^ T 



— — ^1 

C.6S8E CC 

0.1C8E Cl 

“*• I 


4 

4 I 

C.S06E CO 

C.130E Cl 

-1 


4 

4| 

C.SCSE CC 
C.6S8E CC 
C.882E CC 

0.137E Cl 
C.127E Cl 
0.1C9E 01 

^ f 


4 

4 1 

- 1 


4 4 

1 

C.86?E CC 
C.&71E CO 

C.ESeE CC 
c.ecsE CO 
0.636E GO 

— -I 

4 

4 

■ ■ J 

1 

c.esiE CO 

^ 1 


4 4 

I 

C.^3CE CC 

C.72SE CC 

•M. -» ^ 1 . 




C.^85E CC 

0.ES6E 00 




44 I 

C.104E Cl 

0.1C7E Cl 

^ I 



4 4 I 

C-llCE Cl 

C.iiOE Cl 





C.115E Cl 

c,ine Cl 

C.125E Cl 
0.125E 01 




■ -r- - -w M 

4 4 1 

« f 



4 4 1 

C.118E €1 

0.123E 01 


38E 01 


COMFUTEC 

MEASURED 


STRESSES 

STRESSES 


ARE 

ARE 


CIVICEC 0V their peak VAIUE= 5*114H£CTCBAPS 
DIVIDED BY THEIR REAM \ALtE» 8-25CHECTC8AFS 





I 



\ 


0 uo2-ffii * saomHiD 



RUN 9.C5 

tCRFABlSCN CEThttN 
SlRESStS AT .A5*R 
C«C 


CCMPUeO SIRESSeSI^HFRO EKC SmiATlCKSI AKD KEASUREt STRESSES f^HMARCH 11 


PSl 

0 


3C 


€0 


SC 


lic 


15C 


lec 


2K 


2AC 


27C 


30C 


33C 


CCRPtTED 

REASUREC 


I 

-I 
- 1 

-I 

-I 

-1 

-1 

-! 

-I 

-I 

-1 

1 — 

-I 

-I 

1 — 

-I 
- 1 

-1 

-I 

-I 

-I 

-1 

-I 

-1 

-I 

-1 

STRESSES 

STRESSES 


C,ES2E 

1 


CC 


0.136E 01 


, — 

4 

« 

4 

* 




4 


ARE DIVIDED BV IHEIP 
are CIVIOEC BV THEIR 




1 

I 

f^EAh VAtlE* 
l-EAN VAlte^ 


4 4 

♦ 4 

4 4 

A*e04H£CTCB 

?*97CHECTCB 


CGRPUTEC 

MEASURED 

e.ilOE 

Cl 

0.121E 

01 

C.K6E 

Cl 

C<.120E 

Cl 

C.ICAE 

Cl 

C.118E 

Cl 

C*105€ 

Cl 

C.llAE 

01 

C.1C6E 

Cl 

C.1C7E 

01 

0.113E 

Cl 

C.S86E 

CO 

G.llSE 

Cl 

C*SCSE 

CC 

0.122E 

Cl 

C.882E 

CC 

0.U2E 

Cl 

C.S34E 

CC 

C«12CE 

Cl 

C.ICCE 

Cl 

0.115E 

Cl 

C.122E 

Cl 

c.ueE 

Cl 

C.135E 

Cl 

c.ioie 

€1 

C.138E 

Cl 

C.S28E 

CC 

0.127E 

01 

C.858E 

CC 

0.1C3E 

Cl 

C.1S5E 

CC 

C,72^E 

CC 

C-IASE 

CO 

0.A46E 

CO 

C.12AE 

CC 

C,2£SE 

CC 

C.721E 

CC 

C.3C8E 

CC 

0.741E 

00 

o-^sse 

CC 

C.77SE 

CC 

0.7€6E 

CC 

C.62€E 

CC 

0« 1C8E 

Cl 

C.873E 

CO 

0.I28E 

01 

C,SCSE 

CC 

C-13AE 

Cl 

C.S2SE 

CC 

C.12SE 

01 

C,S35E 

CC 

0.1C6E 

Cl 

C-S35E 

CC 

c.ceeE 

CC 

C.S41E 

CC 

C.735E 

oc 

C.Sfc3£ 

CC 

C.713E 

00 

e.lQCE 

Cl 

C.7S5E 

CC 

O-lCtE 

Cl 

C.SACE 

CO 

Q.112E 

Cl 

C.ICSE 

01 

0.117E 

Cl 

C.120E 

Cl 

C.llSE 

Cl 

C,128E 

Cl 

C.118E 

Cl 

0.126E 

01 

C.115E 

Cl 

0.124E 

Cl 


PS 

PS 


TEST5I 



PliK s*ce 

CCKPARISON 6ETMEEN 
STRESSES AT 

C.C 

PSl 1 

30 1 — 

6C 1 — 

SC 1 — 

12C 1 — 

150 --^I — 

180 1 — 

210 1 — 

240 1 — 

27C ! -- 

30C 1 — 

330 1 — 

CCPPUTEO STRESSES 
^^EASUREC STRESSES 


CCMPVTED STRESSESC*MFRCP EPC SIKIATICKSI AAC REASUREC STRESSESI^HNARCH 71 1E$1$> 


C.712E CC 0.1A2E 01 

1 1 
1 4 ♦ 1 

I * ♦ 1 

1 • ♦ ! 

1- ♦ 1 

14-4 1 

1 ♦ 4 I 

— ! 4 4 1 

I + 4 I 

1 + 4 I 

1 4 ♦ I 

I ♦ ♦ 

I 4 4l 

14 ♦ I 

14 4 1 

«^-4-4 1 

4 4 I I 

4 4 1 1 

4-1 1 

4 4 I 

14 I 

1 4 —4 1 

I 4 4 I 

I 4 4 I 

I 44 1 

14 4 1 

4 I 

14 4 ! 

14 4 1 

I r 4 4 — --1 

i 44 1 

I 4*1 

1— 4 4— i 

1 4 4 1 

I 4 4 I 

ARE OlViCEO BT THEIR PEAK VALtE= 4-A48HECTCe ARS 
ARE CIVIOEC 87 THEIR PEAK \ALOE» 7.5SCHECTOE ARS 


CCRPUTEC 
C-S8SE OC 
C.S63E CC 
C.S87E CC 
C.1C3E Cl 
O.IOSE 01 
C,113E Cl 
O.lieE Cl 
C.117E Cl 
C.lllE Cl 
c.nsE Cl 
c.niE Cl 
C.ICCE Cl 
C.S78E CC 
C.E81E CC 
C,717E CC 
c.ee2E CC 
C.C HE CC 
c.eciE CC 
C.eACE CC 
C.726E CC 
C.E3SE CC 
C.S52E CC 
C.ICAE Cl 
C.ICSE Cl 
C.llOE Cl 
C.ICSE Cl 
C. 1C7E Cl 
C.1C7E Cl 
C.ICSE Cl 
C. 112E Cl 
C.115E 01 
C.lllE Cl 
C.llCE Cl 
C.U2E Cl 
C.ICCE Cl 
C.lClE Cl 


MEASURED 
0.124E Cl 
C.118E 01 
O.lllE 01 
0.1C2E Cl 
c.sa^e CC 
C.E35E OC 
C<,7S1E CO 
0.E23E 00 
C.S3SE CC 
O.IUE 01 
G.130E 01 
C.iAlE Cl 
C.140E 01 
0.123E 01 
C.S48E CC 
0.C21E 00 
C.352E CC 
C.228E OC 
erases CO 
CiS15£ OC 
C.830E OC 
C.113E Ci 
G.133E 01 
C.136E €1 
G.128E Cl 
C.llOE Cl 
C-S15E CC 
C.7S4E CC 
C.777E 00 
C.E58E CC 
C.SS8E CC 
C.115E 01 
0.126E Oi 
0.131E 01 
0.132E Cl 
C.12SE Cl 


uo2“Ha s * ooa esroiATHii) 


12. 1C 

CCf'PAPISON e£Tl«££N 
STRESSES AT ,45*R 
C.C 


CCMPITEO STRESSESl^l IFRCR ERC StFtlATlOS) AKC PEASUREC STRESSESC^XHARCH 71 lESlSI 


C.73SE CC 


0 . 


PSI 

c 


3C 


6C 


SC 


12C 


150 


lec 


270 


3CC 


33C 


CCRPUTEC 

REASUREC 


2 1C 


240 


1 

1 — 

-I 
-I 

-I 
-I 

— I — 
-I 
-1 

-I 
•1 

-I 
I 

-I 
-I 

-—I — 

-I 
-I 

I — 

-I 

-! 

-I 

-! 

-I 

-I 

-1 

- I 
-1 

STRESSES 

STRESSES 


4 

* 

♦ 


1 


♦ ♦ 

» 4 

« 4 

4 ♦ 

*• ™ ” ”•■• 4 — A— » — ■ 

4 4 

♦ * 

4 — .—.4 

4 * 

* 4 

4 -. 4 

4 4 ] 

... , 4 4 . .... 

...... 4-*4 

♦ . 

4 , 

— 4 — — . — — — — 

4 
4 

— - 4 .— 4 — 

4 4 

4 4 

♦ . 4 

♦ 4 

-4 * 

4 

4 

-4 ♦ 

4 4 

44 

4 4 
44 

ARE OlVIOEC BY THEIR ^EA^ VAltE» A.604RECTCe ARS 
ARE DIVIDED 8V THEIR ^EA^ VAllE* 7.770HECTCe ARS 


A8E 01 


COMPUTED 
G.llSE 01 
C.llSE Cl 
c.neE Cl 
C.116E Cl 
c.iieE Cl 
C.118E Cl 
c.iieE^ Cl 
C.IISE Cl 
C. IICE Cl 
c.iieE Cl 
0.U3E 01 
C.ICCE Cl 
C.1C2E Cl 
C.S55E CC 
C.B9CE CO 
C.E5CE CC 
C.E2AE CC 
C.EICE CC 
C.E3CE CC 
C.E5CE CC 
C.666E CC 
C.817E CC 
€,t73E CO 
C-eSCE CC 
C.e 32 E 
c.eiiE 
C.EOCE 
C.EOIE GC 
C.835E CC 
c.aeEE 
C.SA3E 
C-IOIE 
C.1C7E 
C.112E 
C.iieE 


CC 

CO 

CC 


CC , 

CC 

Cl 

Cl 

Cl 

Cl 


C.IIEE Cl 


MEASURED 
0.126E Cl 
C.133E Cl 
C.136E Cl 
0.135E Cl 
0.1 2EE Cl 
C.llOE 01 
C.S27E CC 
C.7S3E 00 

CoiesE oc 

C.EEEE CC 
0.1C7E Cl 
C.13CE Cl 
C.IS6E Cl 
C.IASE 01 
0.127E 01 
C.SA7E CC 
C.5S5E CC 
0-331C CC 
C.248E CC 
G.371E CO 
C.CSIE CC 
C-Se3E CC 
C.1256 01 
O.J3eE Cl 
0.12SE Cl 
O.IIOE Cl 
C.EeiE CO 
G.e74E €0 
c.eciE CC 
C.652E CC 
C.7SCE CC 
C.S56E CC 
C.ICSE Cl 
C.IIEE Cl 
0.122E 01 
C.125E Cl 


CD' 


. A 


jPiG ^ 
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UK 12.11 

CKFAPISCN eET«EEN 
iTPtSSES AT ,A 5 *R 
C.C 

PSI 1 

0 1 — 


CCPPneD 51 RESSE Se*HFRCA EPC SIKLATIChS) tht PEASUPEC STRESSES IIWARCH 71 TESTS) 



I 

1 

I 

I >» 

— I*— 

I ♦ 

I 

1 

.^1 

1 > 
! ♦ 

!♦ 

1 * 

.-I — ♦- 

1 

I 


C. 150 E Cl 
I 


XPPUTEC STRESSES 
lEASUREC STRESSES 


I ♦ ♦ 1 

I ♦ ♦ I 

1 *♦ ♦! 

I ♦ ♦ 

I ♦ ♦ 1 

ARE GIVIDEO BV THEIR EEAA VAllE* ^. 592 K£CTCE APS 
ARE DIVIDED 8 Y THEIR PEAK VALUE* 7 . 65 CHECTCE APS 


CCRPUTEC 
C. 135 E Cl 
C. 131 E Cl 
C. 12 AE Cl 
CpIICE Cl 
C.IDTE Cl 
C.ICIE Cl 
C.SEAE CC 
c.seeE CO 

C.ICIE Cl 
C.IC 3 E Cl 
C. 104 E Cl 
C.ICIE Cl 
C.SAAE CO 
C.E 7 CE CO 
C. 802 E CC 
C.TCEE CC 
C. 767 E OC 
C.ECIE CC 
C.S 51 E CC 
C.ES 5 E CC 
C.SIAE CC 
C.SCCE OC 
c.eeiE CO 
C.E 12 E CC 
C. 7 T 5 E CC ' 
C. 76 TE CO 
C. 7 S 3 E CC 
C. 651 E CC 
C. 020 E CC 
C.ICIE Cl 
C.liOE Cl 
C. 11 TE Cl 
C. 12 AE Cl 
C. 129 E Cl 
C. 133 E 01 
C. 135 E Cl 


MEASURED 
0 . 131 E 01 

C. 13 CE Cl 
0 . 123 E 01 

o.iaae ci 

D. 116 E 01 
O.UCE Cl 
C.Sl^E CC 
0 . 7 T 5 E OC 
C. 7 C 5 E CC 
C. 75 CE CO 
C.S 15 E CC 
C.USE Cl 
0 . 135 E 01 
0 . 143 E Cl 
0 . 132 E Cl 
C.ICAE Cl 
o.eTse CO 
0 . 35 CE CC 
C. 1 S 6 E CC 
c.acoE OC 
C. 5 A 3 E CO 
C.S 2 ^E CC 
C. 127 E Cl 
C.H 6 E Cl 
C. 143 E Cl 
0 . 12 GE 01 
C.EE 6 E CC 
C. 61 CE CO 
0.AT6E €0 
C. 53 CE CC 
C. 7 A 3 E CC 
0 .U 3 E Cl 
C. 12 SE Cl 
C. 145 E Cl 
C.ISCE Cl 
C. 14 SE 01 
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12.12 

;Cf*PARlSCM BETktEN CCMPtlEO 


;1BESSE$4*MFRC> ERC SlftLATIChSl AFC FEASUREC STRESSES4-»HRARCH 11 1E$1S> 


STRESSES 

PSl 

C 


30 


AT 

C 


ec 


sc 


120 


150 


180 


210 


2 AC 


27C 


300 


330 


COMPUTED 

PEASUREC 


.A5AR 
C 

1 

-1 

-! 

-1 

-I 

-I 

-1 

1 - 

-I 

-1 

-I 

-I 

— «-I - 
-I 
-1 

1 - 

-1 

-I 

1 - 

-1 

-1 

1 - 

-1 

-1 

1 

-I 

-I 

1 - 

-I 

-1 

“I 

-I 

STRESSES 

STRESSES 


t.EAlE 

1 


cc 


0.16SE 01 


♦ 




1- 

I 
I 

.w-i, 

4 
4 

♦ 1 « 
♦ 

I 

1 

14— 

4 

*I ♦ 

!♦ 

I . 
I~“ 
I 
1 

1 #-. 

4 
1 

1 

♦ 4 
I* 

— I 

1 

♦ I 

I. ^ J ' 




4 

♦ -4- 
4 
4 


4 

"4- 


.♦4- 


ARE 

ARE 


DIVIDED 

DIVIDED 


e\ 

BY 


THEIR 

THEIR 


I 

1 

I 

I 

PEAF 

PEAK 


values 

VALUE* 


4 4 

4 4] 

4 

4 4 I 

44 1 

A-7S2h£CTCEAPS 
7-51CHEC1CEAPS 


COMPUTED 

MEASURED 

C.12SE 

Cl 

C«121E 

Cl 

C.ilEE 

Cl 

c.noE 

Cl 

C.1C5E 

Cl 

0.IC4E 

01 

C.S28E 

CC 

O.SECE 

CC 

C.855E 

tc 

C.SIOE 

00 

C.6A3E 

CO 

0.E31E 

00 

C.CEEE 

cc 

C.153E 

cc 

C.S62E 

cc 

C.733E 

cc 

C.1C3E 

Cl 

C.EICE 

CO 

C.ICEE 

Cl 

e.SElE 

cc 

O.ICAE 

Cl 

O.liOE 

Cl 

t.SCSE 

cc 

C.135E 

Cl 

C.E85E 

CD 

0.138E 

Cl 

C.eiTE 

cc 

C.122E 

01 

C.7S2E 

cc 

C.SCSE 

CC 

C.SISE 

cc 

0.546E 

00 

e.Eeie 

cc 

C.2E4E 

00 

C.SAEE 

cc 

C.170E 

cc 

C.SSIE 

cc 

C.338E 

cc 

c.ssse 

CO 

C.0S4E 

cc 

C.SA2E 

cc 

C.112E 

Cl 

C.E6SE 

cc 

C.HIE 

CI 

C.797E 

cc 

C.ieOE 

01 

C.75EE 

cc 

C.H0E 

Cl 

C.75SE 

cc 

C.llEE 

Cl 

C.EC7E 

cc 

0.i£5E 

oc 

C.EEAE 

cc 

C.A53E 

ec 

C.S7IE 

cc 

0.3416 

cc 

C.IC5E 

Cl 

C.AeiE 

oc 

c.mE 

Cl 

C.174E 

cc 

C.llSE 

Cl 

0.114E 

Cl 

C.125E 

Cl 

C.14SE 

Ci 

C.13CE 

Cl 

C.U7E 

Cl 

C.1356 

01 

C.14CE 

01 

C.138E 

Cl 

C.1S5E 

Cl 

C.I36E 

Cl 

C.126E 

Cl 



% 


. i 
■ f 


T 

f 




• t 


. i 
I 


td U 


s 

• C+- 


I 

I CO 
Ol 


o 



I 

isj 


tt 



Rllfv 12.13 

cc»»p^RiscN eni»£tN 
STRESSES AT .A**R 
O.C 

PS I i 

0 j — 


CCMPITEO S1RESSESI4) IFRCA EPC SI^t^ATIC^S) ARC PEASUREC STRESSES «+l CMARCM 71 TESTS! 


120 


t.?ACE CC 
I 

1—4 

4 * 

44 

«— 4 - 1 — r- 
4 4 I 
44 I 
►4 — 4 — 1 — - — 

4 4 

4 4 


188E 01 




• — 4 — — 


4 1 
44 1 

4 


tCHPUTEC STRESSES 
PEASUREC STRESSES 


4 14 

1 — 4 — - 

4 1 4 

4 4 

— .--♦-I — 

4 1 

* I 

4— 1-*. — 

- 4 4 1 

4 4 

4 14 

4l 4 

1 4 

1 ] 

1 4 

ARE DIVICEC 6> THEIR REAP VAiUE» 
ARE DIVIDED BY THEIR PEAK VAtUE* 


I 

CCPPUTEC 

PEASUREO 

1 

C.126E 

01 

0.1C6E 

01 

I 

C.113E 

Cl 

C.S24E 

CO 


0.^5EE 

CC 

C«E68E 

GO 

J 

C.EITE 

CC 

0.E45E 

CC 

I 

C.7ACE 

CC 

C.SISE 

CO 

! 

0.74AE 

CC 

0.715E 

CC 

1 

c.euE 

CC 

C.742E 

CC 

I 

Co^2AE 

CC 

0.7C7E 

QC 

1 

C.)02E 

Ci 

C.EE5E 

CC 

1 

C.1C7E 

Cl 

C.1C8E 

01 

1 

0.1C5E 

Cl 

0.130E 

01 


C.S76E 

CC 

C.144E 

Cl 

1 

C.893E 

CO 

C.14CE 

Ci 

1 

C.E31E 

CO 

0.117E 

01 

1 

C.EHE 

CC 

C<,EC3E 

CC 

1 

C.E82E 

CC 

0.412E 

CO 

1 

C.S37E 

CC 

C.1S4E 

CO 


C.IOIE 

Cl 

G.141E 

CO 


C.1CA6 

Cl 

0.4C1E 

CC 


C.102E 

Cl 

C.E5CE 

CO 

1 

C.SA3E 

CC 

C.133E 

01 


C.eAAE 

CC 

C-144E 

01 

4 1 

C.759E 

CC 

C.I71E 

Cl 

1 

C*7|EE 

CC 

C.147E 

Cl 

— 1 

C.73tE 

CC 

C.IC3E 

Cl 

I 

C.eOAE 

CC 

0*S5SE 

CO 

I 

C.^C2E 

tc 

C.238E 

CC 

j[ 

C.ICIE 

Cl 

C a i^4E 

CC 


o.noE 

Cl 

0.44CE 

00 

I 

C.U7E 

Cl 

C.EE2E 

CC 

— 1 

C#12AE 

Cl 

C.138E 

01 

4 I 

C.13GE 

Cl 

0.173E 

01 

4 

C.136E 

Cl 

0.1E8E 

Cl 


C . lA IE 

Cl 

c.ieoE 

01 

1 

0.1A2E 

ei 

0.157E 

Cl 


C.138E 

Cl 

G.130E 

01 


A.760HECTCEAFS 

?-540HeCTCEAPS 
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I^UK 14. 1C 
CC|tF/lJ»lSCN 
STRESSES AT 
Q 

PSI 

0 


30 

60 

SC 

12C 

150 

180 

210 

240 

270 

300 


33C 


EETKEEN 

.C 


CCHPLIEO S1RESSES1*MFRC^ EPC SI^tLATIChS) AAC PEASUREC STRESSES 4 «| IttARCH 71 TESTS 


;CRPUTEC STP 
■EASUREC STR 


.773E 

I 


CC 


0. 


«♦ 

>1 ♦ 


« 


ESSES ARE 
ESSES ARE 


CIVIDEC 

DIVIDED 


BT 

BY 


I ♦ * 

1 ♦ 4 

1 4 

I ♦ ♦ 

I ♦ ♦ 

1 ♦ 

14 ♦ 

1---4 

♦ 14 

14 

1-4^ — — 

I 4 

♦ I 4 

!-4 ♦ — — ■ 

14 ♦ . 

1 4 4 

1-4 ♦ 

I 4 ♦ 

1 44 

♦ I-4r-— 

♦ 14 

♦ 14 

— 4,« 

I +4 

I 4 4 

I -——4 

I 4 4 

1 44 

THEIR PEAK VAilE= 5.020HECTCe ARS 
THEIR PEAK VAltE* E.llCHECTGB AES 


55E 01 


CCRPUTE8 


C.117E 

c.neE 

C.115S 

C.114E 

0«114E 

C.115E 

C.116E 

O.liOE. 

C.IHE 

C.119E 

C.1I6E 

0.U2E 

e.lC6E 

C.S36E 

C.EESE 

C.655E 

C.E36E 

C.E31E 

c.ea^E 

C.63SE 
C-E4 2E 
0.E43E 
C.E4CE 
C.E35E 
C.E32E 
C.E34C 
0.E44E 
C-E65E 
C.E97E 
C.S4CE 
C.S^IE 
0.1C4E 
C.1C9E 
0.113E 
C.neE 


01 

u 

Cl 

Cl 

€1 

u 

€1 

Cl 

Cl 

Cl 

Cl 

Cl 

Cl 

CC 

CC 

CC 

CC 

CC 

CC 

CC 

CC 

CC 

CO 

CC 

CC 

00 

CC 

CO 

CC 

CC 

CC 

CC 

Cl 

Cl 

Cl 

Cl 


MEASITREO 
0.121E 01 
C.122E Cl 

c.iasE Cl 
C.1I6E Cl 
C.1S3E 01 
C.U3E 01 
C.^ESE CC 
C.S63E CC 
c.easE CC 

^ C.SIIE- CC 
O.UIE 01 
C.134E €1 
Q«15IE Cl 
0.153E Cl 
C.135S Cl 
0.1C2E Cl 
C.624E CO 
C.3C6E CO 
c.nsE CC 
O.aiTE CC 
C.S72E CC 
C.S49E 00 
C.127E Cl 
0.143E Cl 
o.iacE 01 
C.neE Cl 
C.C76E CC 
o.eaeE co 
0. 528E 00 
C.575E CC 
C.741E CC 
C.S5C6 CC 
C.1J3E Cl 
C.122E 01 
0.125E Cl 
C.123E Cl 



RU^ 14*31 

CC«f/H<ISCN 8nV€£N CCMPllED SIRE SS€ S I «HFRC f £K SIHiL/lKKSI ^i^C RIASUREC- S1BESSES t-^HHARCH' 7 I 1E51S# 
STRESSES AT .a;*R 

CC 


PSl 

0 


3C 


6C 


SC 


C*C 


C.E4CE 

1 

1 

I 

I * 
♦ I 


4 


0«168E 01 
1 

I 

1 

1 

1 


4 14 3 

4 ♦ 1 

1 >» 4 1 

1 ♦ + 

1 4-1 

1 * 4 I 

14 4 1 

1 * 1 

4 * I 

4 4 1 

4 4 1 

* 4 1 

14 4 1 

14 4 1 

4 4 i 

4 1 

-♦ — - — —I 

4 I 

4 4 1 

!>-♦ — I 

4 4 1 

4 4 I 


-1 

12C 1 

-I 
-I 

150 1 

-1 
-1 

lec — I 

-1 
*^ J 

2 10 —I 

-I 
- I 

2AC 1 

-I 

-1 ♦ 

270 1 ♦ 

-I ♦ 

-I 

3CC I 

- 1 
-1 

33C 1 

-1 

-1 

CCRPUTEC STRESSES ARE DIVIDED BY THEIR 


1 

1 
1 

1 

I 

I 

I 

REAR 


MEASURED STRESSES ARE DIVIDED BY THEIR MEAN 


VAtLE* 

VAttE* 


— i 

4 I 

4 4 I 

5*C64HECICBAPS 

7.630HECTCEARS 


COMPUTED 

MEASURED 

C.ilSE 

01 

0.133E 

01 

C.ltCE 

Cl 

C.134E 

01 

C.ICCE 

cx 

0al2S£ 

Cl 

C,S37E 

CD 

C.1I3E 

Cl 

C.S13E 

CC 

C.SC5E 

CC 

C.DSeE 

CC 

c.cecE 

CO 

C.SS3E 

CD 

C*55ie 

00 

C.106E 

Cl 

C.5S2E 

CO 

0. 1I2E 

Cl 

C.S15E 

CC 

0«11SE 

Cl 

0*I15E 

Cl 

0.1I4E 

Cl 

0.1A9E 

01 

C.liOE 

01 

0.IC7E 

Cl 

G.iOAE 

Cl 

0«1C2E 

01 

C*S7SE 

CO 

0«133E 

01 

C.S24E 

CC 

G.CE4E 

oc 

C*66«E 

CC 

C*434E 

GO 

C*fi57E 

CC 

0.13SE 

CO 

C.E3SE 

CC 

C.S55E-CI 

C.82SE 

CC 

C«322E 

CC 

C.E25E 

CC 

0.1 15E 

CC 

C.E31E 

CC 

C.113E 

Cl 

D.846E 

CC 

C.142E 

Cl 

c*E7ce 

CC 

0.I4OE 

01 

C-eSTE 

CC 

C.125E 

Ci 

C.S21E 

CC 

C.iCEE 

Cl 

C*S3€E 

CC 

0.1S2E 

00 

C.S41E 

CC 

c.ccse 

CC 

C.S42E 

CC 

C.5C5E 

CC 

C-SSIE 

CC 

0.1 C5E 

CO 

C.S7SE 

CC 

0.SC3E 

CG 

D.103E 

Cl 

C.llOE 

Cl 

C.llCE 

Cl 

0.U2E 

Cl 

C.I1CE 

Cl 

C.126E 

01 

C.124E 

Cl 

Q.12CE 

01 

0.121E 

Cl 

0.327E 

01 

C.125E 

Cl 

0.I2SE 

Cl 

- 

— 





05 

I ^ CJ 

I 

' i?d 


/r/^ pa , ■ + 
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CCRPARISCN E EThEEN CCHPITEU i IRE SSE S P £K SmiAllChSl AKC PEASUREC STRESSES I^HMARCK 71 lESTSI 

STRESSES AT ,45+R 

-C.2CAE Cl C.C 0.204E 01 

PSI I 1 1 

0 


3C 


6C 


SC 


12C 


ISC 


lec 


210 


240 


270 


3C0 


!3C 


1 

1 



— 1 




-I 

I 

4 4 


-I 

1 

4 4 


— ^ I- 




-! 

I 4 

4 


-I 

I ♦ 

4 


1 -.p.—- 





-1 

1 ♦ 

4 


- 1 

I 

4 4 


1 




•1 

1 

4 

4 

-I 

1 

4 






-I"' 

I 

4 4 


-I 

I 4 

4 






-I 

♦ 1 

4 


-I 

4l 

4 





» « m ^ 

-I 

I 

4 4 


-I 


4 4 





p^*4-.p 

-I 

1 

4 

4 

-I 


4 4 


1 






«I 

I 4 

4 


-1 

I 4 

4 


— - 1 




-I 

I 4 

4 


•I 


4 4 


1- — 




-I 

i 

4 4 


-1 

1 

4 4 


I 





-1 

1 

4 4 


-1 

I 

4 4 



CCRPUTEC STRESSES ARE DIVIOEC BY THEIR PEAK VAllE* 
REASUREC STRESSES ARE DIVIDED BY THEIR PEAK VALUE* 


I 
I 
I 
1 
1 
1 
1 
i 
1 
1 
I 
•# 

1 
I 
1 
1 
I 
1 
I 
1 
I 
I 
I 
I 
1 
I 
1 
1 
I 
] 

I 
1 
I 
1 
1 
1 

5.02CHECTCEAES 

7.750HECTGBARS 


COMPUTED 

REASURED 

C.llCE 

Cl 

G.13CE 

Cl 

C.S2CE 

CC 

C.123E 

01 

c.ieee 

cc 

0*1€3E 

Ci 

C.t€5E 

CC 

c.Tses 

cc 

0.E64E 

cc 

0.419E 

00 

C,7S5E 

cc 

C.237E 

cc 

(.SICE 

cc 

C.313E 

cc 

C.106E 

Cl 

C.571E 

CO 

C.123E 

Cl 

C.IC3E 

Cl 

C.I31E 

ci 

C.154E 

Cl 

C.132E 

Cl 

C.1?4E 

Cl 

C.12EE 

Cl 

C.2C4E 

Cl 

C.IJEE 

Cl 

o.ncE 

Cl 

ColCBE 

Cl 

C.124E 

0! 

Cossee 

cc 

c.sesE 

cc 

C.^35E 

cc 

-0.578E-C2 

C«€€7E 

cc 

-C.274E 

cc 

C.E4SE 

cc 

-CpISCE 

cc 

CpEUE 

cc 

C.31SE 

cc 

C.798E 

cc 

0.547E 

cc 

C.7S4E 

cc 

C.152E 

Cl 

CoElSE 

cc 

C.1E4E 

Cl 

C.652E 

cc 

G.I82E 

01 

C.SC2E 

cc 

C.ISCE 

Cl 

C.S47E 

cc 

0.IC3E 

Cl 

C.S74E 

cc 

0.SS3E 

cc 

C.SIEE 

cc 

C.343E 

cc 

c.sice 

cc 

0.343E 

cc 

C.S67E 

cc 

C.557E 

cc 

C.SSCE 

cc 

' C.BCIE 

cc 

C.103E 

€1 

C.12CE 

Cl 

G.114E 

Cl 

C.H2E 

01 

C.124E 

Cl 

C.151E 

Cl 

C.131E 

Cl 

C.1S2E 

Cl 

0.132E 

Cl 

O.HBE 

01 

C.124E 

u 

C.H3E 

Cl 


CP 

-4 


I 

t?3 
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CCWf^RISCN eEl^EEN CCKPtTEt S IRE SS£ S ii IFRCF EPC StM;L«nCRS} 4 RC Pt/iS<i 5 <£e SIRESSE 3 
STRESSES AT , 4 ‘*R 

- 0 .E 2 SE Cl t.C 0 . 229 E 01 


'»I 4 »ARCH T 1 1 ESTSI 


♦ — 

* 

♦ 

4. 

♦ ♦ 

♦ ♦ 



-1 ♦ I 4 I 

-I ♦ I * 1 

18 C 1 !-♦ ’ 1 

-1 I 4 I 

-1 I ♦ ♦ I 

210 1 “’■*£ ♦ “““♦ J 

-1 1441 

- ! I * * * 

24 C I I 4-4 I 

-1 I ♦ 4 . I 

-! 1 ♦ 4 1 

27 C I * !“♦ -- 4 — — 1 

-I 144 1 

- ! 14 4 1 

3 CC I — --I- 44 : ♦! 

-1 1441 

- I I 4 4 1 

■3‘3Q 1 — — — — -» 4 — - 4 — 1 

" -I I 441 

- 1 1 44 1 

CCRPUTED STRESSES ARE DIVIDED BV IHEIP REA^ \ALtE= 5 . 112 KECTCEAFS 
WEASUREC STRESSES ARE DIVIDED BT THEIR REAk VALL'E= T.T 50 HECTC 6 ARS 


CCRRUTED 

HEASUREO 

C.USE 

Cl 

0 . 1 35 E 

01 

C.S 1 SE 

cc 

C. 114 E 

01 

C-ISEE 

CO 

C.E 73 E 

CO 

C.SE 2 E 

cc 

C. 205 E 

CC 

C. 5 C 8 E 

cc 

0 . 36 SE 

CC 

C. 553 E 

cc 

C. 12 CE 

CO 

C. 7 C 2 E 

cc 

C. 23 SE 

cc 

C.S 15 E 

cc 

0 .S 87 E 

00 

C. 112 E 

(1 

C. 114 E 

Cl 

c.iasE 

Cl 

0 . 173 E 

Cl 

C. 138 E 

Cl 

0 . 217 E 

01 

c.nei 

Cl 

0 . 22 CE 

Ci 

C. 123 E 

01 

C. 2 C 0 E 

01 

0 . 123 E 

Cl 

0 . 136 E 

01 

C. 113 E 

Ci 

C. 0 C 7 E 

CC 

C. 104 E 

Cl 

- 0 . 5 E 4 E-C 1 

C.S 55 E 

cc 
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C.S 83 E 
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C.^TSE 
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C. 154 E 
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Cl 

C.UCE 

01 
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Cl 

C.U 5 E 

Cl 

C. 139 E 

Cl 

c.ieoE 

Cl 

C. 134 E 

Cl 

C. 15 CE 

Ci 
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RUN 16.08 
COMPARISOVJ 
STRESSES A 


BETWEEN COMPUTED STRESSES C ♦)( FROM EMC SIMULATIONS! AND MEASURED STRESSES! ♦MMARCH 71 TESTS! 
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COMPUTED 
0.107E 01 
0.107E 01 
0. 106E 01 
0.109E 01 
0.112E 01 
0.116E 01 
0.I21E 01 
0.125E 01 
0.127E 01 
0.126E 01 
0.123E 01 
O.110E 01 
O.lllE 01 
0.104E 01 
0.979E 00 
0.934E 00 
0.907E 00 
0.896E 00 
0.895E 00 
0.899E 00 
0.902E 00 
0.899E 00 
0.888E 00 
0.871E 00 
0.850E 00 
0.828E 00 
0.812E 00 
0.805E 00 
0.812E 00 
0. 834E 00 
0.871E 00 
0.917E 00 
0.966E 00 
O.IOIE 01 
0.104E 01 
0.106E 01 


MEASURED 
0.1I3E 01 
0.123E 01 
0.128E 01 
0.124E 01 
0.112E 01 
0.980E 00 
0*901E 00 
0.937E 00 
0.109E 01 
0.131E 01 
O.l^BE 01 
0.133E 01 
0.139E 01 
0.108E 01 
0.TO5E 00 
0.383E 00 
0.226E 00 
0.288E 00 
0.536E 00 
0«871E 00 
0.117E 01 
0.132E 01 
0.129E 01 
O.lllE 01 
0.B95E 00 
0.708E 00 
0.651E 00 
0i727E 00 
0.887E 00 
0.105E 01 
0.116E 01 
0.117E 01 
O.lllE 01 
0.103E 01 
0.997E 00 
0.104E 01 
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RUN 16.10 

COMPARISQ'4 BETWE6N CO'^PUTEO STRESSES! ♦!* FROM EMC SIMULATIONS! AND 
STRESSES AT .45#R 
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COMPUTED STRESSES ARE DIVIDED BY THEIR MEAN VALUE= 4.620HECTOBARS 


MEASURED STRESSESt^HMARCH 71 TESTS! 


COMPUTED 
O.lllE 01 
O.lllE 01 
0.1 13E 01 
0.117E 01 
0.120E 01 
0.120E 01 
0.117E 01 
0.1 lOE 01 
0.102E 01 
0.937E 00 
0.890E 00 
0.893E 00 
0.946E 00 
0.103E 01 
0.I12E 01 
0.118E 01 
0. 120E 01 
0.1186 01 
0.113E 01 
0.1066 01 
O.IOOE 01 
0.953E 00 
0.91 IE 00 
0.868E 00 
0 . 81 ^ 00 
0.7606 00 
0.709E 00 
0.6816 00 
0.690E 00 
0.742E 00 
0.8296 00 
0.930E 00 
0.102E 01 
0.109E 01 
O.lllE 01 
0.1126 01 


MEASURED 
0.117E 01 
0.12^6 01 
0.130E 01 
0.130E 01 
0.I23E Oi 
0.1126 01 
0.995E 00 
0.915E 00 
0.9096 00 
0o978E 00 
0.1096 01 
0.119E 01 
0.123E 01 
0.1186 01 
0.1056 01 
0.895E 00 
0.7676 00 
0.7176 00 
0.7616 00 
0.875E 00 
O.IOIE 01 
0.109E 01 
O.llOE 01 
0;l02E 01 
0.882E 00 
0-7^76 00 
0.666E 00 
0.663E 00 
0.7376 00 
0.852E 00 
0.9656 00 
O-IO^^E 01 
0.1076 01 
0ol07E 01 
0.1086 01 
O.lllE 01 
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